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PREFACE 

This  report  describes  the  research  completed  as  a  result  of  the 
contract  F49620-81-K-0D10  awarded  to  Beth  Israel  Hospital,  Boston, 
Massachusetts  by  the  Air  Force  Office  of  Scientific  Research  (AFOSR). 
The  contract  technical  monitor  (CTM)  during  the  development  of 
experimental  equipment  was  Col.  Irving,  and  during  the  later  part  of 
the  project  was  Dr.  Berry. 

The  delay  in  the  completion  of  this  contract  was  due  to  two 
principal  reasons: 

1)  The  development  and  construction  of  the  multi-axial  materials 
test  system  required  for  the  testing  of  spine  segments  was  more 
complex  and  time  consuming  than  anticipated. 

2)  Administrative  problems  placed  constraints  on  the  budget  which 
delayed  hardware  development  and  limited  personnel. 

An  Interim  report  was  submitted  to  the  contract  technical  monitor 
in  August,  1981  describing  the  design  analysis  of  test  hardware. 

The  authors  acknowledge  the  contribution  of  several  individuals 
who  participated  in  this  work.  The  assistance  of  Joe  Lelli  is 
acknowledged  for  his  help  during  the  early  development  of  test  system 
hardware.  Y.  K.  Kou  is  gratefully  acknowledged  for  design  of  the 
central  interface  between  the  testing  system  and  the  computer. 
Further  appreciation  is  due  Dr.  Berry  for  his  understanding  and 
patience  while  this  report  was  being  prepared. 

The  support  provided  by  the  National  Institutes  of  Health  through 
Grant  AM33066,  and  by  a  Grant  from  General  Motors  Research  Laboratory 
after  the  expiration  of  this  contract  made  possible  the  completion  of 
the  work  reported  here. 
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1.  INTRODUCTION 


1.1  BACKGROUND 

Reduction  of  the  risk  of  injury  to  the  cervical  region  of  the 
spine  resulting  from  extreme  loading  conditions  of  the  type  observed 
in  motor  vehicle  accidents  and  emergency  ejection  from  military 
aircraft  is  of  concern  to  many.  Detailed  data  of  the  relationship 
between  loads  and  displacements,  the  vertebral  strength,  and  the  modes 
of  injury  in  the  spine  are  required  to  define  the  tolerance  of  the 
spine  to  injury. 

Experimental  methods  previously  used  for  in  vitro  biomechanical 
studies  of  the  spine  have  been  limited.  Methods  developed  to  simulate 
the  complex  combination  of  loads  and  displacements  in  the  spine 
generally  employ  small  loads  in  comparison  to  the  normal  physiologic 
and  fracture  conditions.  Other  methods  for  study  of  spinal  fracture 
employ  uniaxial  material  testing  machines.  A  primary  goal  of  this 
project  was  to  develop  a  multiaxial  materials  testing  system  capable 
of  simulating  the  load  displacement  conditions  for  sagittal  motions 
and  load  specimens  to  failure.  This  new  multiaxial  testing  system  was 
used  to  investigate  the  mechanics  of  the  cervical  spine  and  the 
tolerance  of  the  cervical  spine  to  injury. 


1.2  PROGRAM  SUMMARY 

The  objectives  of  this  research  were:  1)  to  construct  a  multiaxial 
spinal  testing  apparatus  for  measuring  the  material  properties  of  the 
spine,  2)  to  determine  the  geometry  of  human  cervical  spine  segments, 
3)  to  measure  the  viscoelastic  load-displacement  properties  of  the 
cervical  spine,  and  4)  to  model  the  cervical  spine.  A  multiaxial 
testing  apparatus  was  constructed  with  three  independent  components  of 
displacement  input  and  six  independent  components  of  load  measurement. 
This  system,  called  the  Planar  Testing  Apparatus  (PTA),  can  generate 
planar  motions  needed  to  characterize  the  sagittal  response  of  spine 
segments.  A  study  of  the  viscoelastic  properties  of  two  vertebrae 
lumbar  spine  segments  was  conducted  following  the  completion  of  the 
PTA  to  check  and  demonstrate  the  test  system.  The  results  from  seven 
lumbar  specimens  were  also  included  in  this  report. 

Seven  cervical  spine  segments  were  tested  both  nondestructively 
and  destructively  to  determine  the  relationships  between  loads  and 
displacement,  vertebral  strength,  and  modes  of  failure.  The  Planar 
Testing  Apparatus  was  used  to  measure  the  mechanical  properties  of 
specimens  by  applying  controlled  loads  or  displacements  in  the 
sagittal  plane.  The  spine  segments  consisted  of  three  cervical 
vertebrae  and  their  interconnecting  soft  tissue,  discs,  and  ligaments. 
Three  specimens  were  taken  from  the  upper  cervical  region,  one  from 
the  central  (C3-C5)  region,  and  three  from  the  lower  cervical  spine. 
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A  noninvasive  electro-mechanical  displacement  measuring  apparatus  was 
developed  which  monitored  the  six  degree-of-f reedom  motion  of  the 
middle  vertebral  body  as  it  moved  relative  to  the  two  adjacent  fixed 
bodies.  Kinematic  equations  were  derived  to  compute  the  three- 
dimensional  translations  and  the  Eulerian  rotations  at  the  center  of 
the  moving  vertebra.  The  results  included  in  this  report  are 
presented  in  three  parts:  A)  the  development  of  the  Planar  Testing 
Apparatus,  B)  the  viscoelastic  properties  of  lumbar  spine  segments  in 
compression  (a  pilot  study),  and  C)  the  viscoelastic  and  strength 
properties  of  cervical  spine  segments. 


2.  PLANAR  TESTING  APPARATUS  DEVELOPMENT 

The  Planar  Testing  Apparatus  (PTA)  is  a  multi-degree-of-f reedom 
material  testing  machine  (Fig.  1).  The  concept  and  kinematic  design 
were  originated  by  W.T.  Edwards  and  Y.F.  Kou  at  the  Orthopaedic 
Biomechanics  Laboratory  and  developed  under  this  grant.  The  interface 
of  the  testing  machine  to  the  PDP  11/03  microcomputer  and  the  control 
protocol  of  the  testing  machine  were  also  developed  during  this 
project . 

The  general  layout  of  the  Planar  Testing  Apparatus  is  similar  to 
most  electrohydraulic  material  testing  machines  (such  as,  INSTRON 
model  1331)  in  that  it  contains  a  moving  test  stage  and  a  stationary 
load  cell.  The  test  specimens  installed  on  the  test  stage  can  be 
loaded  independently  in  the  vertical  direction  (tension/compression, 
z-axis),  the  horizontal  direction  (shear,  x-axis),  and  rotated  about 
the  axis  perpendicular  to  the  vertical  and  horizontal  axis  (y-axis). 
Figure  2  shows  the  orientation  of  the  coordinate  axes  with  respect  to 
the  test  stage.  For  this  project,  functional  spinal  unit  specimens 
were  placed  on  the  test  stage  such  that  x-z  plane  represented  the 
sagittal  plane  and  positive  x  was  in  the  anterior  direction  of  the 
specimen  (Fig.  3) 

2.1.  HARDWARE 

2.1.1.  Load  Frame 

Test  specimens  are  mounted  on  the  Planar  Testing  Apparatus  load 
frame  between  a  moving  test  stage  below  and  a  stationary  load  cell 
above.  The  load  cell  is  bolted  to  a  76x127x940  mm  cross  head  which  is 
mounted  on  two  61  mm  diameter  stainless  steel  poles.  The  i_^osshead  is 
clamped  onto  each  of  the  stainless  steel  poles  by  2  sets  of  nuts  and 
bolts.  The  crosshead  can  slide  freely  on  the  poles  when  the  nuts  and 
bolts  are  loosened,  but  is  secured  by  tightening  the  nuts  and  bolts. 
Specimens  up  to  914  mm  in  length  and  254  mm  in  width  can  be  mounted  on 
the  Planar  Testing  Apparatus  for  testing.  The  maximum  dimensions  of 
the  Planar  Testing  Apparatus  are  2000  mm  in  height,  660  mm  in  width 
and  952  mm  in  length. 
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2.1.2.  Load  Cell 

A  multi-component  load  cell  (Robert  A.  Denton,  Inc.  Model  1687 
Serial#  75)  was  mounted  onto  the  crosshead.  This  load  cell  is  a 
cylinder  of  140  mm  in  diameter  and  76  mm  in  height.  The  load  cell  can 
resolve  a  maximum  of  17,792  N  in  the  z  direction  (tension- 
compression),  a  maximum  of  8,896  N  in  the  x  direction  (anterior- 
posterior  shear)  and  a  maximum  of  1,129  N-m  about  the  y  axis  (sagittal 
bending).  Signal  conditioning  for  the  load  cell  is  obtained  with 
Vishay  2120  strain  gage  conditioners. 


2.1.3.  Test  Stage 

The  test  stage  (Figure  4)  is  a  140x254  mm  moving  platform 
supported  from  below  by  three  linear  hydraulic  actuators.  The  three 
hydraulic  actuators  are  arranged  so  that  two  of  the  actuators  are 
attached  to  one  axle  of  the  test  stage  and  the  third  actuator  attached 
to  the  remaining  axle.  Tapered  roller  bearings  are  used  on  the  axles 
to  reduce  any  play  in  the  actuators  and  test  stage  attachment 
mechanism.  Extreme  care  was  taken  during  the  assembly  of  the  test 
stage  and  bearings  to  properly  align  the  three  actuators  in  the  plane 
of  motion.  Holes  are  drilled  and  taped  on  the  surface  of  the  test 
stage  for  installation  of  test  fixtures  to  support  the  specimens. 

2.1.4.  Hydraulic  Actuators 

The  hydraulic  actuators  used  on  the  Planar  Testing  Apparatus  have 
three  main  components. 

(a)  the  servovalve; 

(b)  the  piston;  and 

(c)  the  LVDT  position  sensor. 

The  servovalve  responds  to  the  incoming  electrical  current  and 
directs  the  oil  flow  path  to  extend  or  to  retract  the  actuator  piston. 
The  position  of  the  hydraulic  actuator  is  sensed  by  a  linear  variable 
differential  tranducer  (LVDT)  which  is  installed  within  the  actuator's 
housing.  A  M00G  power  oscillator  (M00G  123-105)  is  used  to  provide 

excitation  for  the  LVDT. 

Three  high  performance  M00G  (M00G,  Inc.,  East  Aurora,  N.Y.) 
cylinder-servovalve2  assemblies  (model  A085-233)  were  used.  This 
actuator  has  710  mni  piston  crossectional  area  and  152  mm  stroke.  The 
servovalve  is  rated  at  207  kPa  maximum  supply  pressure  with  a  maximum 
no-load  velocity  of  1244  mm  per  second  (at  137  kPa  supply)  and  an 
equivalent  first-order  time  constant  of  0.002  second  (M00G  Inc.  :  M00G 
A085  Servoactuators ,  Data  sheet  no.  850-382). 


2.1.5.  Electro-hydraulic  Controller 

MOOG  122-104  analog  controllers  are  used  to  control  the  position 
of  the  hydraulic  actuator.  A  demodulator  is  built  into  the  controller 
to  condition  the  A.C.  signal  from  the  LVDT,  providing  D.C.  signal 
proportional  to  the  actuator  position.  The  D.C.  signal  indicates  the 
position  of  the  hydraulic  actuator  and  serves  as  the  feedback  in  the 
control  circuit.  The  inputs  to  the  closed-loop  control  circuit  are 
the  mean  level  voltage  and  the  command  voltage.  The  mean  level 
voltage  offsets  the  positions  of  the  actuator  so  that  the  full  range 
of  the  actuator  can  be  used  (Figure  5).  Once  the  mean  level  voltage 
has  been  set,  it  remains  unchanged  during  the  test  sequence.  The 
command  voltage  is  the  voltage  which  actually  controls  the  motions  of 
the  test  stage.  The  control  block  diagram  of  the  MOOG  122-104  analog 
servo  controller  is  shown  in  Figure  6. 

2.1.6.  Display  Panel 

A  display  panel  allows  manual  control  of  the  Planar  Testing 
Apparatus,  (Figure  7).  The  panel  contains  three  identical  sets  of 
knobs,  indicators  and  circuits  which  correspond  to  the  three  hydraulic 
actuators.  Four  LEDs  are  used  to  indicate  the  range  setting 
(100%, 50%, 10%, 5%) .  The  100%  range  corresponds  to  a  maximum  actuator 
range  of  motion  of  +75  mm.  Digital  voltmeters  (DVM)  are  used  to 
display  the  position  of  the  hydraulic  actuators  in  percent  of  the 
maximum  displacement  of  the  selected  range.  Selection  switches  are 
used  to  select  the  DVM  display  options.  The  display  on  the  digital 
voltmeter  can  be  the  input  command  signal,  the  mean  level  setting,  the 
LVDT  feedback  signal,  the  error  signal  or  the  actuator's  limit 
setting.  This  information  is  also  displayed  in  percent  of  the  maximum 
actuator  displacement  of  the  selected  range. 

In  addition  to  displaying  the  status  of  the  hydraulic  actuator, 
the  display  panel  also  interprets  the  command  signal  and  the  range 
selection  signal  from  the  POP  11/03  microcomputer.  The  range  selection 
signal  activates  the  analog  switches  in  the  circuit  and  directs  the 
command  signal  to  appropriate  potentiometers  where  the  incoming  signal 
is  dropped  to  the  proper  range.  The  resulting  signal  from  the  panel's 
input  to  the  MOOG  controller.  The  upper  and  lower  limits  of  the 
actuator  displacement  can  be  set  on  the  panel.  When  a  pre-selected 
limit  is  exceeded,  the  reed  relay  on  the  MOOG  D.C.  level  detector 
(Model  123-103)  shuts  off  the  hydraulic  supply  to  the  actuators.  The 
mean  level  voltage  is  also  set  on  the  display  panel.  Figure  8  shows 
the  signal  flow  diagram  of  the  Planar  Testing  Apparatus. 

2.1.7.  PDP  11/03  Microcomputer 

A  PDP  11/03  microcomputer  (Digital  Equipment  Corp.,  Maynard,  MA) 
is  used  both  as  a  calculating  and  a  controlling  device.  The  functions 
of  the  microcomputer  used  by  the  Planar  Testing  Apparatus  including 
start  up,  range  selection,  range  calibration,  open-loop  test  stage 
position  control,  closed-loop  load  control  and  data  acquisition.  The 
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microcomputer  is  set  up  to  take  18  channels  of  data.  The  maximum  data 
rate  is  10  data  points  per  second  for  displacement  control  and  7  data 
points  per  second  for  load  control.  The  control  cycle  rate  is  100 
Hertz  and  30  Hertz  for  displacement  control  and  load  control 
respectively.  A  block  diagram  of  the  complete  system  is  shown  in 
Figure  9. 

2.2.  GEOMETRIC  CONSIDERATION 

Since  the  Planar  Testing  Apparatus  is  capable  of  three  degree-of- 
freedom  motions,  a  reference  point  on  the  x-z  plane  (sagittal  plane) 
is  necessary  to  define  the  position  and  the  orientation  of  the  test 
stage.  Figure  10  shows  a  simplified  drawing  of  the  test  stage  and  the 
hydraulic  actuators.  Points  E  and  D  are  the  two  axles  of  the  test 
stage  and  Points  A,  B  and  C  are  the  three  axles  on  the  base  frame  to 
which  the  three  hydraulic  actuators  are  attached.  Parameters  ES  and 
TS  are  input  by  the  user  to  define  the  reference  point  T.  The 
reference  point  T  is  the  center  of  rotation  and  represents  the  test 
stage  position  on  the  x-z  plane.  The  origin  0  is  defined  at  the  mid¬ 
point  between  A  and  C.  The  position  of  point  T  on  the  x-z  plane 
(with  respect  to  the  origin  0),  and  the  angle  theta  (defined  by  DE  and 
the  horizontal  line),  define  the  spatial  position  of  the  test  stage. 
This  configuration  is  referred  to  as  the  absolute  orientation. 

2.3.  SYSTEM  PRECISION 

The  main  factor  affecting  the  precision  of  the  Planar  Testing 
Apparatus  system  is  the  resolution  of  the  A/D  and  D/A  modules  on  the 
PDP  11/03  microcomputer.  The  PDP  11/03  is  a  16-  bit  machine  and  the 
MINC  software  uses  the  4  most  significant  bits  as  status  bits,  leaving 
12  bits  to  resolve  the  analog  voltage  signal.  The  full  range  of  the 
D/A  module  is  +  10V  with  a  resolution  of  0.0048V  per  bit  The  full  range 
of  the  A/D  module  is  +5. 1175V  with  a  resolution  of  0.0025V  per  bit. 
The  M00G  122-104  servo-controllers  have  been  modified  such  that  the 
LVDT  reading  of  +5V  indicates  the  full  stroke  (+76.2  mm)  of  the 
hydraulic  actuators  and  the  command  voltage  of  +10V  moves  the 
hydraulic  actuator  +75  mm  at  100X  range  setting.  This  arrangement 
allows  the  highest  possible  resolution  with  the  available  equipment. 
The  PDP  11/03  can  thus  detect  a  change  in  actuator  displacement  of 
0.037  mm  and  can  control  an  actuator  displacement  of  0.0018  mm  at  5£ 
range  setting  with  all  actuators  near  mid-stroke. 

Assuming  that  the  time  constant  of  the  electronic  control 
circuitary  is  negligible  compared  to  the  time  constant  of  the 
hydraulic  actuator  and  that  the  inertia  of  the  actuator  is  negligible, 
the  system  can  be  represented  by  the  first  order  system  shown  in 
Figure  11.  The  gain  of  the  servocontroller  (K)  can  be  adjusted  from 
between  1  to  20  with  the  present  circuit  design.  From  experiment,  the 
rise  time  of  the  hydraulic  actuator  with  no  mass  attached  is  0.06 
second  with  a  unit  step  voltage  input.  The  system  is  stable 
throughout  the  complete  range  of  the  gain  adjustment.  With  the  test 
stage  mounted  on  the  hydraulic  actuators,  the  added  mass  of  the  test 
stage  cannot  be  neglected  and  the  system  becomes  second  order.  Hence, 


the  stability  is  affected  by  the  gain  setting  [24].  Different  gain 
settings  have  been  tried  and  gain  settings  under  10  were  found  to  be 
stable. 

2.4.  RANGE  OF  MOTION  AND  LOAD  CAPACITY 

The  configuration  of  the  Planar  Testing  Apparatus  has  been 
optimized  to  meet  the  motion  range  and  load  capacity  required  for 
tests  in  spinal  mechanics  [10].  The  kinematic  design  minimizes  the 
gross  actuators  motion.  The  maximum  horizontal  range  is  330  mm,  the 
maximum  horizontal  force  is  36,474  N,  the  maximum  vertical  force  is 
42,256  N  and  the  maximum  moment  is  31,610  N-m.  The  ranges  of  motion 
and  the  load  capacity  of  the  test  stage  at  different  orientation  are 
tabulated  in  Tables  1  to  3  respectively. 

2.5.  SOFTWARE 

2.5.1.  Organization 

The  control  program  is  organized  on  three  levels.  The  main 
program  contains  the  menu  of  all  the  available  options  for  the  control 
of  the  apparatus.  The  second  level  is  the  control  programs.  All  the 
command  signals  are  generated  at  this  level.  Data  acquisition  and 
storage  also  take  place  at  this  level.  The  third  level  includes  the 
utility  programs.  This  level  manages  interfacing  between  the  software 
and  the  laboratory  hardware. 

2.5.2.  Start  Up  Routine 

A  range-calibration  provides  the  voltage  relationship  between  the 
D/A  output  voltage  and  the  voltage  actually  seen  by  the  hydraulic 
controller.  There  are  four  range-calibration  data  files  which 
correspond  to  the  four  ranges  (100%, 50%, 10%, 5%)  used  for  hydraulic 
control.  During  start  up,  one  of  the  four  range  calibration  data 
files  is  read  into  the  computer  program.  At  the  same  time,  a  data 
file  defining  actuator  location  vs  the  LVDT  feedback  voltage  is  also 
read  in.  A  voltage  is  then  output  to  the  hydraulic  controller  so  that 
the  error  signal  becomes  zero  (Figure  12).  This  voltage  is  simply 
governed  by  the  operational-amplifier  equation  : 


0.0=(E  /R  )  +  (E  /R  )  +  (E.,,nT/R.,.nT) 

com  com  mean  mean  LVDT  LVDT 

where, 

E  :  voltage 
R  :  resistance 
com  :  command 
mean  :  mean  level 
LVDT  :  LVDT  feedback 


The  appropriate  output  voltage  (E  )  can  be  calculated  from 
Equation  (1)  and  output  through  the  D/S°module  to  the  MOOG  hydraulic 
controller,  after  the  LVDT  feedback  voltage  (ELVnT^  anc*  t*ie  mean  level 
voltage  (E  )  are  read  into  the  computer.  Error  indicating  LEDs  on 
the  displayefanel  go  dark  when  it  is  safe  to  actuate  the  hydraulic 
actuator . 

2.5.3.  Open-Loop  Control  Routine 

The  role  of  the  computer  in  the  control  scheme  defines  open  or 
closed  loop  control.  An  open-loop  control  routine  is  used  to  position 
the  test  stage.  The  user  can  input  the  desired  translational  and 
angular  displacement  or  input  the  final  absolute  orientation  of  the 
test  stage.  To  facilitate  testing,  portions  of  this  information  can 
also  be  input  through  a  data  file.  A  ramp  function  is  generated  by 
the  control  program  to  move  the  test  stage  from  the  current 
orientation  to  the  specified  orientation  and  the  corresponding  output 
voltage  stored  in  an  integer  array.  A  test  sequence  of  ten  different 
orientations  can  be  specified  and  the  routine  generates  functions  to 
command  the  test  stage  to  the  positions  in  sequence.  The  voltage 
information  array  is  reused  for  data  storage  to  reduce  the  core  memory 
required  to  run  the  control  routine.  Nevertheless,  a  rate  of  10  data 
samples  per  second  can  still  be  attained  for  18  channels  of  data.  The 
open-loop  routine  can  also  use  the  recorded  data  file  to  reproduce  the 
displacement  history  from  the  particular  test. 


2.5.4.  Closed-Loop  Control  Routine 

The  closed-loop  control  routine  changes  the  orientation  of  the 
test  stage  so  that  the  load  (tension-compression,  shear  and  bending 
moment)  specified  by  the  user  can  be  obtained.  This  routine  uses  an 
integral  control  scheme  (Figure  13)  so  that  small  disturbances  do  not 
affect  the  final  load  value.  The  control  algorithm  makes  a  uni¬ 
directional  search  until  the  loads  reach  the  pre-determined  values. 
The  toe  region  of  the  load-displacement  response  for  biological  tissue 
does  not  affect  the  machine  in  searching  for  the  loads.  The  feedback 
path  is  completed  by  the  three-component  load  cell  and  the  signal 
conditioner  (Figure  8).  The  user  is  required  to  supply  the  desired 
loads,  the  test  stage  velocity  limit  and  the  feedback  gain.  The 
magnitude  of  the  feedback  gain  is  dependent  upon  the  stiffness  of  the 
specimen.  For  higher  stiffnesses,  a  lower  feedback  gain  is  required 
to  maintain  stability.  Since  the  routine  uses  an  integal  control 
scheme,  a  high  gain  can  result  in  load  oscillation  [24J. 

A  linearized  version  of  the  test  stage  position  control 
calculation  is  used  to  reduce  computation  time  during  each  closed-loop 
cycle.  The  sampling  frequency  (sweeping  rate)  of  the  closed-loop 
control  cycle  is  set  at  30  Hz.  The  control  strategy  employed  assumes 
no  coupling  effect  between  different  components  of  loads. 
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2.5.5.  Calibration  Routine 

There  are  four  potentiometers  in  each  of  the  three  sets  of 
circuits  in  the  display  panel.  These  potentiometers  reduce  the  D/A 
signal  to  one  of  the  four  range  settings  (100%, 50%, 10%, 5%)  selected  by 
the  user.  The  calibration  routine  serves  to  calibrate  the 
potentiometer  settings  to  ensure  an  accurate  relationship  between  the 
D/A  voltage  output  and  the  voltage  actually  seen  by  the  hydraulic 
controller.  The  routine  utilizes  the  existing  range  calibration  data 
to  move  the  hydraulic  actuators  to  different  points  along  the  whole 
displacement  range  of  the  actuator.  The  LVDT  feedback  voltages  at 
these  points  are  recorded.  By  using  Equation  (1),  the  voltage  seen  by 
the  hydraulic  controller  can  be  calculated.  A  linear  regression  is 
then  used  to  relate  the  D/A  module  output  to  the  voltage  calculated. 
The  result  is  then  stored  in  a  new  range-calibration  data  file. 

2.5.6.  Set  Mean  Level  Routine 

The  mean  level  voltage  is  used  to  offset  the  position  of  the 
individual  hydraulic  actuators  so  that  an  experiment  can  be  performed 
at  the  most  convenient  region  of  the  stroke  range.  Figure  5  shows  a 
simplified  illustration  of  the  operation  of  the  mean  level  voltage. 
The  displacement  shown  in  Figure  5  can  be  translational  or  rotational 
displacement.  Since  the  Planar  Testing  Apparatus  has  three  degrees  of 
freedom,  the  mean  level  setting  for  the  three  hydraulic  actuators 
cannot  be  randomly  selected.  This  routine  allows  the  users  to  set  the 
mean  level  voltage  such  that  a  zero  command  voltage  corresponds  to  the 
current  orientation  or  to  any  specified  absolute  orientation  of  the 
test  stage. 

2.5.7.  Test  Stage  Position  Adjustment  Routine 

Adjustment  of  the  test  stage  position  is  necessary  during 
installation  of  the  specimen.  The  adjustment  routine  independently 
moves  the  test  stage  in  the  z,  the  x  directions  or  rotates  the  test 
stage  about  the  y  axis.  The  adjustment  has  fine  and  coarse  control 
and  can  be  terminated  instantly  by  the  users  during  the  execution  of 
the  routine. 

2.6.  SUMMARY  OF  PLANAR  TESTING  APPARATUS 

The  Planar  Testing  Apparatus  provides  a  new  approach  to 
biological  material  testing.  In  the  past,  tests  of  multi-axial, 
coupled  spinal  mechanics  generally  used  static  weights  applied  with 
cables  and  pulleys.  The  displacement  was  measured  by  mechanical  or 
electro-mechanical  devices.  For  equipment  of  this  type,  there  is  no 
effective  way  to  control  the  displacements  of  the  specimen  in  more 
than  one  axis. 


The  ability  of  the  Planar  Testing  Apparatus  to  perform  load  or 
displacement  control  in  a  mul ti-degree-of-f reedom  environment  not  only 
makes  more  realistic  physiological  testing  possible,  but  has  opened 
several  new  areas  of  interest.  For  example,  the  response  of  the 
intervertebral  joint  to  different  combinations  of  displacement 
directions  and  rates  can  be  evaluated.  It  should  now  be  possible  to 
evaluate  the  stiffness  matrix  of  the  intervertebral  joint  directly 
instead  of  inverting  the  flexibility  matrix  [25].  With  the  use  of  the 
multi-component  load  cell,  the  coupling  effect  can  also  be  evaluated. 


3.  VISCOELASTIC  PROPERTIES  OF  THE  LUMBAR  SPINE 

Following  the  development  of  the  Planar  Testing  Apparatus,  a 
brief  project  was  undertaken  to  refine  the  testing  capabilities  of 
this  system.  The  lumbar  region  of  the  spine  was  selected  for  this 
study  because  the  specimens  were  available  for  testing  and  the 
cervical  spine  specimens  could  be  in  this  way  saved  for  the  main  part 
of  the  study,  presented  in  the  next  section  of  this  report,  and 
because  the  biomechanical  properties  of  the  lumbar  spine  are  relevant 
to  the  better  understanding  of  both  low  back  pain  and  injuries  to  the 
spine. 


3.1  REVIEW  OF  THE  LITERATURE 

Compressive  loads  applied  to  the  spine  are  supported  anteriorly 
by  the  intervertebral  disc  and  the  longitudinal  ligaments,  and 
posteriorly  by  the  facet  joints  and  the  posterior  ligaments.  The 
intervertebral  disc  and  the  facet  joints  in  the  lumbar  region  of  the 
spine  have  both  been  identified  as  possible  sites  of  mechanically 
induced  pain.  Excessive  mechanical  stress  can  produce  pain  and  injury 
in  the  anterior  or  the  posterior  spine  components.  Increased  anterior 
load  can  promote  disc  degeneration,  disc  herniation,  or  disc  proplase. 
Increased  posterior  loading  may  produce  soft  tissue  injuries  or 
osteoarthritis  of  the  facet  joints.  An  understanding  of  the  load 
distribution  characteristic  of  the  intervertebral  joint  is  thus  an 
important  step  in  understanding  the  etiology  of  low  back  pain  and  the 
mechanisms  of  injury  in  the  spine. 

Low  back  pain  is  one  of  the  most  frequent  and  costly  medical 
complaints  in  this  country  [13,33,34].  Kelsey,  et  al.  [13]  showed 
that  the  back  and  spine  disorders  were  the  most  frequent  causes  of 
activity  limitation  in  persons  under  age  45  with  a  decrease  of  28.6 
days  a  year  in  work  capacity  per  100  subjects  age  from  25  to  44  in  the 
United  States.  Low  back  pain  affects  more  than  half  of  the  working 
population  at  some  time  during  their  working  careers  [33].  Based  on 
data  from  237  long-term  male  employees  just  before  their  retirement  at 
age  62-  65,  Rowe  found  that  56X  of  workers  experience  low  back  pain 
sufficient  to  require  medical  treatment  [30].  Drivers,  material 
handlers  and  office  workers  are  amoung  the  highest  risk  groups  for  low 
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back  pain  [33,34].  Even  though  drivers  and  office  workers  do  not 
usually  participate  in  heavy  lifting,  both  kinds  of  jobs  involve  long 
period  of  sitting.  Therefore,  the  time  dependent  load  distribution 
characteristic  of  the  intervertebral  joint  cannot  be  neglected  when 
considering  the  cause  of  low  back  pain. 

In-vivo  measurements  indicate  that  large  forces  occur  in  the 
lumbar  spine  during  the  activities  of  daily  living.  Nachemson 
measured  interdiscal  pressure  in-vivo  and  compared  the  result  with  the 
pressure  obtained  from  cadaver  specimens  [22].  He  pointed  out  that 
the  pressure  in  the  L3-L4  disc  corresponds  to  a  force  as  high  as  2.5 
times  body  weight  when  the  subject  is  sitting  or  standing  at  20 
degrees  flexion.  His  experiment  also  established  that  the  interdiscal 
pressure  is  related  to  the  mechanical  stress  applied  to  the 
intervertebral  joint  [23].  Since  the  load  supported  by  the  lumbar 
disc  is  high,  it  has  been  hypothesized  that  mechanical  stresses  play 
an  important  role  in  the  pathogenesis  of  low  back  pain.  However,  the 
specific  mechanisms  involved  are  poorly  understood  [37]. 

Both  the  anterior  and  posterior  portions  of  the  intervertebral 
joint  display  nonlinear,  coupled  viscoelastic  load-displacement 
properties  [29].  We  may  hypothesize  that  a  difference  in  the 
viscoelastic  properties  of  the  disc  and  the  facets  can  cause  a  re¬ 
distribution  of  load  between  these  functional  spinal  unit  components 
in  response  to  sustained  compressive  load.  Thus,  when  a  constant  axial 
compressive  load  is  applied  to  the  intervertebral  joint,  the  relative 
load  distribution  between  the  intervertebral  disc  and  the  facet  joints 
vary  with  time. 

The  objectives  of  this  investigation  of  the  lumbar  spine  were  to 
assess  the  time  dependent  load  distribution  characteristics  of  the 
intervertebral  joint  using  physiological  combinations  of  displacements 
and  loads.  The  specific  objectives  were: 

1)  Measure  the  creep  response  of  intact  lumbar 
functional  spinal  units; 

2)  Determine  the  fraction  of  total  applied  load 
supported  by  the  disc  and  the  facets;  and 

3)  to  generate  an  idealized  viscoelastic  model  for  the 
lumbar  functional  spinal  units. 

The  advanced  testing  methodologies  needed  for  these  biomechanical 
measurements  were  developed  for  the  Planar  Testing  Apparatus  and 
demonstrated  using  this  system. 


18 


3.2  MATERIALS  AND  METHODS 

The  load-displacement  and  creep  characteristics  of  7  lumbar 
functional  spinal  units  (FSU)  were  studied.  Physiologic  load 
conditions  were  selected  to  test  the  distribution  of  compressive  load 
between  anterior  and  posterior  functional  spinal  unit  components  over 
time. 

3.2.1.  Pretest  Preparation 

Two  vertebrae,  lumbar  functional  spinal  units  (FSU)  of  L1-L2  and 
L3-L4  were  used  in  this  experiment.  Human  spines  were  obtained  from 
donated  cadavers  from  the  Harvard  Medical  School  Anatomical  Gifts 
Program.  These  specimens  were  collected  within  48  hours  after  death 
and  were  placed  into  tightly  sealed  plastic  bags  to  prevent 
dehydration  and  were  frozen  at  -20  degrees  Celsius  until  needed.  X- 
rays  were  taken  of  each  specimen  to  screen  out  those  with  metastatic 
bone  disease,  previous  spinal  fractures  or  other  disorders  known  to 
alter  the  mechanical  properties  of  the  spine.  Prior  to  testing,  the 
spines  were  thawed  over  night  at  room  temperature.  Excess  soft 
tissues  were  cleaned  from  the  specimens  leaving  the  interconnecting 
ligaments  intact.  The  spines  were  then  divided  into  L1-L2  and  L3-L4 
functional  spinal  unit  specimens  by  dissection  of  the  adjacent  discs. 
The  exposed  discs  at  the  T12-L1,  L2-L3  and  L4-L5  levels  were  graded 
from  0  (normal)  to  3  (degenerated)  using  the  grading  method  of  Galante 
191. 

Seven  lumbar  specimens  were  used  for  this  experiment.  The  age, 
sex,  disc  grade,  and  causes  of  death  for  these  cadaver  specimens  are 
shown  in  Table  4a. 

3.2.2.  Test  Fixture  and  Mounting  Protocol 

A  special  test  fixture  was  designed  for  the  load-  displacement 
and  creep  tests  Figure  14.  With  this  fixture,  a  constant  reference 
configuration  was  exactly  reproducible  during  repeated  removal  and 
installation  of  the  fixture  in  the  Planar  Testing  Apparatus.  The  two 
identical  parts  of  the  fixture  were  constructed  from  a  rectangular,  25 
mm  thick  aluminium  plate.  An  aluminium  ring  31.8  mm  in  height  with  an 
outer  diameter  of  133.3  mm  and  12.7  mm  wall  thickness  was  attached  to 
the  center  of  each  plate.  Four  63.5  mm  long  set  screws  positioned  90 
degree  apart  were  screwed  through  the  wall  into  the  aluminium  ring 
radially  from  the  outside.  Another  four  socket  head  screws  were 
screwed  into  the  wall  along  the  axial  direction  of  the  ring  to  secure 


the  aluminium  ring  onto  the  plate. 


Epoxy  was  used  to  seal  the  gap 


between  each  ring  and  plate  so  that  there  was  no  leakage  through  the 
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superior  and  inferior  spinous  processes  of  each  vertebra.  Half  of  the 
inferior  vertebral  body  of  the  specimen  was  then  placed  into  the  cup 
on  the  lower  plate  of  the  fixture.  By  using  the  long  set  screws 
around  the  aluminium  cup,  the  specimen  was  carefully  positioned  so 
that  the  mid-plane  of  the  intervertebral  disc  was  parallel  to  the 
lower  plate.  Resin  ('Castolite'  Buehler  Ltd.  No.  20-8120-009)  was 
then  poured  into  the  aluminium  cup  until  the  resin  covered  the  long 
metal  screw  in  the  inferior  vertebral  body,  the  wood  screw  in  the 
spinous  process,  and  the  four  set  screws  of  the  ring,  all  of  which 
served  to  secure  the  spinal  specimen  within  the  resin.  The  resin 
hardened  in  about  1.5  hours. 

The  specimen  was  kept  moist  at  all  times.  For  mounting,  wet 
cotton  padding  was  wrapped  around  the  specimen  and  periodically 
injected  with  water  by  a  syringe  to  maintain  a  moist  environment. 
During  the  test,  an  additional  plastic  saranite  film  was  placed  around 
the  specimen.  Radiographs  were  used  to  verify  proper  positioning  of 
the  spinal  specimens  in  the  test  fixture,  Figure  15.  Four  threaded 
rods  were  screwed  superior-inferiorly  into  the  four  corners  of  the 
lower  plate  and  four  aluminium  tubes  of  equal  length  (82.5  mm)  were 
sleeved  over  each  rod.  The  upper  plate  was  then  placed  onto  the 
fixture  with  the  cup  facing  down  toward  the  specimen  and  the  four 
threaded  rods  through  the  four  corresponding  holes  on  the  upper  plate. 
The  upper  plate  was  supported  by  the  four  equal  length  aluminium  tubes 
and  became  parallel  to  the  lower  plate.  Nuts  were  put  onto  the  rods 
to  secure  the  two  plates  in  position.  The  length  of  the  aluminium 
tubes  were  selected  so  that  about  half  of  the  superior  vertebral  body 
was  embedded  in  the  cup  on  the  upper  plate.  The  superior  vertebral 
body  was  held  in  place  by  the  set  screws  of  the  upper  aluminium  cup 
while  resin  was  poured  into  the  cup  to  secure  the  superior  vertebral 
body. 


After  the  resin  was  set,  the  lower  and  upper 
the  threaded  rods  and  aluminium  tubes  were  moun 
and  the  load  cell  of  the  Planar  Testing  Appara 
functional  spinal  unit  was  positioned  on  the  tes 
Testing  Apparatus  so  that  the  sagittal  plane  o 
the  same  plane  as  the  plane  of  motion  (x-z) 
Apparatus.  This  initial  orientation  of  the  test 
as  the  set  up  orientation. 

3.2.3.  Neutral  Point 


plates  together  with 
ted  to  the  test  stage 
tus  (Figure  16).  The 
t  stage  of  the  Planar 
f  the  specimen  was  in 
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Most  biological  tissues  (including  the  intervertebral  joint)  are 
viscoelastic  and  non-linear.  Figure  17  shows  a  typical  load- 
displacement  response  of  a  functional  spinal  unit  for  one  complete 
cycle  of  displacement.  The  flat  region  around  the  zero  load  point  has 
been  called  the  toe  region  or  the  neutral  zone.  Several  different 
definitions  have  been  given  for  the  neutral  zone  and  the  neutral  point 
(8,29).  In  our  tests,  the  neutral  point  for  each  functional  spinal 
unit  was  defined  as  the  mid-point  between  the  two  intersections  of  the 
hysteresis  curve  with  the  zero  load  line  (points  A  and  B  in  Figure 
17).  The  neutral  point  and  neutral  zone  were  defined  for  the  three 
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components  of  sagittal  displacement  by  measuring  the  load-displacement 
response  for  each  independent  motion.  The  neutral  orientation  was 
defined  as  the  position  of  the  test  stage  where  all  three  components 
of  displacement  (axial  translation  :  x,  anterior-posterior  translation 
:  z,  and  sagittal  rotation  :  Qy)  were  located  at  their  neutral  points. 
The  initial  translation  (x=0,  z=0)  and  initial  rotation  (theta=0)  for 
all  experiments  on  each  specimen  were  also  defined  at  the  neutral 
orientation. 

After  the  specimen  was  placed  in  the  Planar  Testing  Apparatus  and 
the  set  up  orientation  defined,  complete  cycles  of  motions  in  x,z 
translation  and  sagittal  rotation  were  executed  on  the  testing  machine 
to  generate  hysteresis  curves  for  the  specimen.  The  three  load- 
displacement  (tension-compression,  translational  shear  and  sagittal 
bending)  curves  were  plotted  on  the  CRT  of  the  computer.  The 
displacements  from  the  set  up  orientation  to  the  neutral  orientation 
were  estimated  from  the  CRT  display  and  recorded.  The  test  stage  was 
then  moved  to  this  neutral  orientation  and  the  orientation  stored  into 
memory  for  the  remainder  of  the  tests. 

3.2.4.  Functional  Spinal  Unit  Geometry 

For  these  experiments,  the  superior  and  inferior  vertebral  bodies 
were  assumed  to  be  rigid.  The  dimensions  of  the  set  up  of  the 
specimen  and  fixture  were  obtained  from  the  lateral  radiograph  (Figure 
15).  The  magnification  factor  of  the  radiograph  was  calculated  by 
comparing  the  actual  dimensions  between  two  markers  placed  on  the 
fixture  and  its  images  on  the  X-ray.  The  conventions  used  to  define 
the  geometry  of  the  superior  and  inferior  vertebral  bodies  are  shown 
in  Figure  18.  The  center  of  the  vertebral  body  was  defined  as  the 
intersection  of  the  two  lines  connecting  the  mid-points  of  AB,  CD  and 
AD,  BC.  The  loads  applied  to  the  specimen  were  resolved  to  the  center 
of  the  superior  vertebral  body.  The  center  of  rotation  of  the  test 
stage  was  placed  at  the  center  of  the  inferior  vertebral  body. 

3.2.5.  Experimental  Protocol 

The  experiment  was  divided  into  five  parts.  Part  I  involved 
defining  the  neutral  point  and  setting  up  the  experiment.  The  load- 
displacement  responses  of  the  intact  functional  spinal  unit  and  the 
isolated  anterior  components  were  measured  in  Parts  II  and  IV, 
respectively.  The  creep  responses  of  the  intact  function  spinal  unit 
and  the  anterior  portion  of  the  function  spinal  unit  were  measured  in 
Parts  III  and  V. 

The  geometry  of  the  functional  spinal  units  was  measured  from  the 
radiographs,  and  the  data  is  presented  in  Table  4b. 

Part  I  -  Neutral  Point  Location 

After  the  fixture  with  the  intact  spinal  specimen  was  installed 
onto  the  test  stage,  the  aluminium  tubes  and  the  threaded  rods  were 
removed.  The  x,z  coordinates  and  the  sagittal  angle  were  recorded  as 
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the  set  up  orientation.  Starting  from  the  set  up  orientation,  the 
specimens  were  cycled  through  compression  and  tension,  anterior  and 
posterior  shear  (at  a  rate  of  0.5  mm  per  second),  and  sagittal  bending 
(at  a  rate  of  1  degree  per  second).  Each  component  of  motion  was 
applied  separately  so  that  three  series  of  measurements  were  made. 
The  translations  and  rotation  between  the  set  up  orientation  and  the 
neutral  orientation  were  determined  from  the  load-displacement  curves. 
The  test  stage  was  then  moved  to  the  neutral  orientation. 

Part  II  -  Hysteresis  Curves 

Starting  from  the  neutral  orientation,  the  specimens  were 
subjected  to  compression  and  tension,  and  anterior  and  posterior  shear 
at  a  rate  of  0.5  mm  per  second  and  to  sagittal  bending  at  a  rate  of  1 
degree  per  second.  These  tests  generated  hysteresis  curves  in  the  z 
,x  directions  and  sagittal  rotation  respectively.  The  hysteresis 
curves  were  used  to  verify  the  material  properties  and  served  as  a 
reference  for  other  parts  of  the  experiment. 


Part  III  -  Intact  Functional  Spinal  Unit  Testing 

This  part  of  the  test  used  the  Planar  Testing  Apparatus  to 
simulate  a  constant  axial  load  applied  to  the  intact  functional  spinal 
unit.  A  mixed  load-displacement  control  option  was  used.  The 
functional  spinal  unit  specimens  were  loaded  in  compression  while 
controlling  the  anterior-posterior  shear  load  at  zero  and  the  flexion- 
extension  angle  at  a  constant  setting.  The  final  axial  compressive 
load  was  set  at  1200  N  using  a  displacement  ramp  function  limited  to 
0.1  mm  per  second.  The  anterior-  posterior  shear  load  was  set  at  0  N 
and  the  rate  of  sagittal  rotation  was  set  at  0  degrees  per  second.  A 
feedback  gain  of  0.001  mm  per  N  was  chosen  for  the  experiments.  For 
each  specimen,  tests  were  performed  at  fixed  flexion-extension  angles 
of  0  degrees,  2  degrees  extension  and  2  degrees  flexion  in  sequence. 
The  test  lasted  1  hour  for  each  setting.  After  each  test,  the  test 
stage  was  returned  to  the  neutral  orientation  and  the  specimen  was 
relaxed  for  1  hour  between  tests.  The  displacement  history  throughout 
this  part  of  the  experiment  was  recorded  by  the  computer  and  stored 
into  data  files  for  future  use. 

Following  these  tests  of  intact  functional  spinal  units,  the  four 
aluminium  tubes  and  threaded  rods  were  re-installed  and  the  fixture 
together  with  the  specimen  were  removed  from  the  test  stage  of  the 
Planar  Testing  Apparatus.  The  intact  specimen  was  then  stored  in  the 
cold  room  at  four  degrees  Celsius  for  12  hours  until  the  next  part  of 
the  experiment. 

Part  IV  -  Isolated  Intervertebral  Disc 

The  load-displacement  curves  of  the  intervertebral  disc  were 
measured  in  this  part  of  the  test.  The  specimens  were  reinstalled 
into  the  Planar  Testing  Apparatus.  The  aluminium  tubes  and  threaded 
rods  were  removed  from  the  fixture  and  the  test  stage  moved  to  the 


neutral  orientation  obtained  from  Part  I.  The  load-displacement  in 
tension-compression,  translational  shear  and  sagittal  bending  (Part 
II)  were  repeated  and  the  data  compared  with  previous  data.  The  test 
fixture  was  then  removed  from  the  test  sfage  and  the  resin  casting 
removed  from  the  aluminium  fixture  by  unscrewing  the  set  screws  around 
the  aluminium  ring.  The  posterior  elements  with  the  embedding  material 
were  then  transected  by  cutting  vertically  through  the  pedicles.  The 
intervertebral  discs  were  then  placed  back  into  the  fixture  and  the 
set  screws  were  used  to  obtain  the  previous  fixture  dimensions. 
Automobile  body  putty  was  then  used  to  fill  the  gap  left  open  by  the 
removed  embedding  material.  The  fixture  together  with  intervertebral 
disc,  the  aluminium  tubes  and  threaded  rods  were  re-installed  onto  the 
testing  machine.  The  test  stage  was  moved  to  the  neutral  orientation 
based  on  the  data  obtained  previously.  Hysteresis  curves  were 
generated  in  x,  z  directions  and  in  sagittal  rotation  about  the  y 


Part  V  -  Creep  History  Test  -  Isolated  Anterior  Components 

After  the  response  of  the  isolated  anterior  components  was 
obtained  (Part  IV),  the  displacement  history  recorded  in  Part  III  was 
used  to  deform  the  isolated  intervertebral  disc.  The  load  resulting 
from  this  deformation  was  recorded  and  stored  for  later  comparisons 
with  the  intact  specimens.  The  testing  sequence  of  0  degrees  flexion- 
extension,  2  degrees  extension  and  2  degree  flexion,  and  the 
relaxation  time  of  one  hour  between  tests  were  the  same  as  in  Part 
III. 


3.2.6.  Analysis 

The  data  transducers  allowed  simultaneous  recording  of  the  three 
components  of  load,  the  positions  of  the  three  hydraulic  pistons  and 
time.  Six  channels  on  the  PDP  11/03  computer  were  used.  A  data 
acquisition  rate  of  0.25  seconds  per  data  point  was  used  during  the 
load  transient  and  a  rate  of  30  seconds  per  data  point  was  used  during 
creep.  After  each  test,  the  data  obtained  were  stored  on  floppy 
diskettes.  The  files  were  then  transferred  from  the  diskettes  to  the 
VAX  11/750  computer  for  further  analysis.  From  these  six  channels  of 
data,  translations  in  the  x,z  directions  and  the  sagittal  plane 
rotation  as  well  as  the  three  components  of  load  were  calculated. 
These  data  were  then  read  into  the  RS/1  statistics  package  for 
plotting  and  function  fitting  (RSI  Bolt,  Beranek  and  Newman  Research 
Systems,  Cambridge,  Massachusetts). 

The  x,z  displacements  and  the  sagittal  angle  were  defined  as  zero 
at  the  neutral  orientation.  The  three  components  of  load  reading  at 
the  neutral  orientation  were  used  as  an  offset  and  were  substracted 
from  the  corresponding  component  of  load  reading  so  that  the  loads 
always  read  zero  at  the  neutral  orientation. 


Equation  (2)  was  used  to  represented  the  response  of  two  parallel 
3-parameter  Kelvin  solid  models  (Fig.  19).  Our  experiment  used  a  ramp 
displacement  of  0.1  mm/sec  until  the  axial  compressive  load  of  1200  N 
was  obtained.  Hence,  a  standard  step  input  response  function  could 
not  be  used  to  fit  the  experimental  data.  In  this  case,  the  data  was 
fit  from  the  initial  time  (t  when  the  constant  axial  load  (1200  N) 
was  just  obtained  and  the  term  representing  the  initial  conditions  was 
incorporated  into  the  step  response  function  equation.  From  the 
equation,  it  is  clear  that  the  effect  of  the  initial  condition 
gradually  diminishes.  In  fact,  the  values  of  the  constants  (A,B,C,D) 
calculated  using  the  step  response  function  were  very  close  to  the 
values  of  the  constants  calculated  from  Equation  (2). 

S=A( l-exp( -( t-to)/B)+C(l-exp(-(t-to)/D)+E(exp(-(t-tQ)/B)+ 
(So-E)exp(-(t-to)/D)  (2) 

where , 

A,B,C,D,E  =  constant 

S  =  displacement 

t  =  time 

t  =  initial  time 

S  =  initial  displacement 

o  v 

For  our  test  conditions,  the  response  described  by  this  equation  is 
equivalent  to  the  response  of  a  system  with  two  times  constants  and  an 
initial  elastic  response,  a  5-parameter  solid.  Because  of  this 
similarity,  this  model  will  be  referred  to  as  the  5-parameter  model 
below.  Equation  2  was  fit  to  the  displacement  versus  time  curve  of 
the  intact  intervertebral  joint  using  the  Marquardt-Levenberg  method 
from  the  RSI  statistics  package.  This  technique  combines  the 
advantages  of  the  high  convergence  rate  of  the  Gauss-Newton  procedure 
and  the  simple  alogrithm  of  the  method  of  steepest  descent  in 
calculation  of  the  least-square  solution. 

When  the  displacement  history  of  the  intact  intervertebral  joint 
(Part  III)  was  reapplied  to  the  transected  specimen  (Part  V),  the 
recorded  load  represented  the  portion  of  load  supported  by  the 
intervertebral  disc  during  the  previous  test  with  the  intact 
intervertebral  joint.  The  difference  between  the  constant  applied 
load  and  the  load  recorded  from  the  transected  specimen  was  then  taken 
as  the  load  supported  by  the  facet.  Since  the  data  acquisition  rate 
for  these  two  experiments  ( I I I , V )  was  the  same.  The  data  could  be 
processed  point  by  point,  i.e.  without  interpolation. 


3.3.  RESULTS 


The  results  of  this  study  of  the  lumbar  FSU  show  that  the  axial 
compression  tests  at  0  degrees  flexion-extension  angle  were  most 
consistent  throughout  the  experiment.  The  data  presented  here  focuses 
on  the  results  for  the  0  degree  flexion  flexion-extension  angle, 
although  data  for  the  2  degrees  flexion  and  2  degrees  extension  tests 
are  also  shown.  The  results  of  the  load  distribution  experiments  are 
presented  here  in  the  following  order: 

(a)  the  hysteresis  characteristic  of  the  functional 
spinal  units, 

(b)  the  displacement  versus  time  for  the  intact 
functional  spinal  unit; 

(c)  the  axial  load  across  the  intervertebral  disc 
versus  time; 

(d)  the  percent  of  total  applied  load  supported  by  the 
disc  versus  time;  and 

(e)  the  load-displacement  of  the  intervertebral  disc 
during  the  load  transient. 


3.3.1.  Hysteresis 

Figure  20.1-5  shows  the  hysteresis  curve  of  the  functional  spinal 
unit  specimens.  These  were  obtained  before  the  first  series  of  load 
control  creep  tests  (Part  III)  and  just  before  the  displacement 
history  test  (Part  V).  The  hysteresis  curves  are  used  to  compare  the 
material  properties  of  the  functional  spinal  units. 

3.3.2.  Displacement  versus  Time  -  Intact  Functional  Spinal 
Unit 

Functional  spinal  unit  specimens  (at  0  degrees  flexion- 
extension)  under  a  constant  axial  load  of  1200  N  displayed  the  typical 
behavior  of  a  Kelvin  solid.  Seventy  percent  of  the  total  displacement 
observed  over  the  hour  of  testing  occurred  within  the  first  30  seconds 
of  the  test.  This  initial  response  was  followed  by  a  slower  creep 
period  which  continued  throughout  the  remainder  of  the  loading  period. 
The  displacement  versus  time  plot  of  the  intact  functional  spinal 
units  are  shown  in  Figures  21.1-5.  The  measured  data  of  the  0  degrees 
flexion-extension  tests  were  fit  with  the  Kelvin  solid  model  Equation 
(2),  using  the  Marquardt-Levenberg  method.  These  best  fit  curves  are 
shown  in  Figures  22.1-5  along  with  the  experimental  data. 

There  are  two  time  constants  in  the  5-parameter  Kelvin  solid 
model,  Equation  2.  For  7  specimens,  the  long  time  constant  had  a  mean 
of  2620.0  seconds  (S.D.  769  seconds)  and  the  short  time  constant  had 
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a  mean  of  79.6  seconds  (S.D.  _  32.6  seconds).  The  results  of  the  curve 
fitting  are  tabulated  in  Table  5.  For  the  five-parameter  Kelvin  solid 
model,  the  data  were  fit  from  the  initial  time,  t  where  the  axial 
force  just  attained  a  steady  load  at  1200  N.  The  initial  conditions 
of  the  displacement  at  time,  t  were  incorporated  into  the  fitting 
procedure. 

3.3.3.  Intervertebral  Disc  Load 

The  displacement  during  the  closed-loop  cycle  in  Part  III  of  the 
experiment  was  compared  with  the  open-loop  cycle  in  Part  V  to  confirm 
that  the  displacement  of  the  intact  functional  spinal  unit  was 
repeatable  in  Part  V  of  the  test,  The  displacements  were  reproduced  to 
within  0.05mm.  Typical  results  are  plotted  in  Figure  23. 

Figures  24.1-5  show  the  axial  load  supported  by  the 
intervertebral  disc  at  all  three  sagittal  rotation  angles  by 
replicating  the  motions  of  the  intact  specimen  from  Part  III.  The  disc 
load  level  of  the  2  degrees  flexion  test  is  consistently  higher  than 
that  of  the  2  degrees  extension  test.  However,  the  disc  load  level  of 
these  two  tests  were  not  consistently  higher  or  lower  than  the  0 
degrees  flexion-extension  tests.  The  disc  load  level  of  the  0  degrees 
flexion-extension  tests  varied  from  about  1000  N  to  1150  N  where  as 
the  2  degrees  flexion  tests  and  the  2  degrees  extension  tests  varied 
from  400  N  to  1400  N. 

Figures  25.1-4  show  the  time  dependent  characteristics  of  the 
axial  load  across  the  isolated  intervertebral  disc  at  0  degrees 
flexion-extension  angle.  The  load  on  the  anterior  components  is 
expressed  as  a  percentage  of  the  total  applied  axial  load.  The 
transient  response  during  the  first  30  seconds  of  the  test  was  not 
consistent  between  specimens.  The  load  supported  by  the 
intervertebral  disc  alone  may  range  from  50%  to  150%  during  the 
transient  period  following  the  application  of  load.  Thus,  the  first 
30  seconds  of  data  is  not  shown  in  Figure  25.  The  remaining  data  were 
fit  by  a  linear  least  squares  method  and  are  also  shown  in  the 
figures.  For  4  specimens,  the  mean  initial  loads  supported  anteriorly 
was  77.2%  (S.D.  9.7%)  of  the  total  intact  load.  The  mean  load  after 
one  hour  was  88.8%  (S.D.  _  4.0%)  nf  the  total  intact  load  (specimens 
B,D,G,,  Table  6.  The  gradual  increase  in  the  anterior  load  was 
observed  for  all  but  one  (specimen  F)  functional  spinal  unit  specimen. 
Stiffening  of  the  disc  was  observed  with  specimen  C.  In  this  case, 
the  load  supported  by  the  intervertebral  disc  was  higher  than  the 
applied  load.  This  violates  static  equilibrium.  Hence,  specimen  F 
and  C  were  not  used  in  the  calculation  of  the  mean  load  and  its 
standard  deviation. 


26 


3.4  DISCUSSION 

3.4.1.  Planar  Testing  Apparatus 

In  previous  studies,  application  of  unidirectional  loads  were 
always  accompanied  by  motions  in  more  than  one  direction  due  to 
coupling  effect.  Panjabi,  et  al.  [25,26]  measured  the  flexibility 
matrix  of  thoracic  functional  spinal  units  and  noted  a  significant 
coupling  motion  between  axial  and  anterior-posterior  translation. 
Schultz,  et  al.  found  that  sagittal  bending  and  anterior-posterior 
shear  are  coupled  in  the  lumbar  spine  [31].  Lin,  et  al.  [17]  and  Yang 
and  King  [38]  used  wedge  shape  knife  edge  to  apply  eccentric  load  to 
spinal  specimens.  Adams  and  Hutton  [1]  Cyron  and  Hutton  [5]  use 
rollers  and  plates  to  simulate  flexion  and  extension.  These  load 
application  methods  inevitably  generate  anterior-posterior  shear 
during  axial  compression.  The  Planar  Testing  Apparatus  developed  for 
this  project  was  able  to  eliminate  coupling  forces  and  generate  a 
unidirectional  load  of  physiological  magnitude  (1200  N).  The 
advantages  of  the  Planar  Testing  Apparatus  over  other  testing  systems 
were  demonstrated. 

In  these  lumbar  spine  experiments,  a  mixed  load  and  displacement 
control  scheme  was  used  to  control  anterior  and  posterior  shear.  The 
center  of  rotation  was  defined  at  the  center  of  the  inferior  vertebral 
body.  When  the  spinal  specimen  is  extended  or  flexed,  anterior- 
posterior  shear  is  generated  due  to  the  deformation  of  the 
intervertebral  disc.  In  extension,  a  posteriorly  directed  shear 
force  is  exerted  on  the  inferior  vertebral  body.  In  the  same  fashion, 
specimens  under  flexion  generate  an  opposite  response.  The  Planar 
Testing  Apparatus  maintained  a  zero  anterior-posterior  shear  by 
translating  the  inferior  vertebral  body  posteriorly  during  extension 
tests  and  by  translating  anteriorly  during  flexion.  The  typical 
anterior-posterior  translation  (x)  of  the  specimen  needed  to  produce 
zero  anterior-posterior  loads  under  flexion  and  extension  conditions 
is  shown  in  Figure  26.  In  this  was  a  pure  unidirectional  compressive 
load  was  produced. 

The  displacement  resolution  of  the  testing  machine  can  affect  the 
stability  of  the  closed-loop  load  control  of  the  Planar  Testing 
Apparatus.  At  100%  range  setting,  the  resolution  of  the  hydraulic 
actuator  is  0.036  mm.  When  there  is  an  error  in  load  detected,  a 
minimum  change  in  actuator  length  of  0.036  mm  is  necessary.  The 
actual  translation  in  the  x,  z  directions  and  sagittal  plane  rotation 
varies  depending  on  the  test  stage  position.  For  stiff  materials, 
small  changes  in  displacement  produce  large  changes  in  load,  and  thus 
may  result  in  load  oscillation.  The  resolution  of  the  linear 
displacement  and  the  rotational  displacement  of  the  Planar  Testing 
Apparatus  at  10%  range  setting  (0.0036  mm  per  bit)  proved  to  be  stable 
in  our  experiments. 

The  gain  can  also  affect  the  stability  of  load  control,  Figure 
11.  High  gain  can  induce  load  oscillation  and  a  low  feedback  gain 


increases  the  rise  time  of  the  system  [24].  The  gain  setting  of  0.001 
was  selected  for  stability  and  reasonable  rise  time  (approximately  12 
seconds)  for  the  one  hour  experiments. 

One  control  cycle  is  defined  as  the  process  to  acquire  data,  to 
calculate  and  to  execute  the  command.  A  high  control  cycle  frequency 
is  required  to  avoid  large  step  change  in  displacement.  Large  step 
change  in  displacement  can  cause  "overshoot"  in  the  motions  of  the 
test  stage.  Even  though  the  three  components  of  load  are  coupled, 
this  coupling  is  reasonably  small,  and  the  control  algorithm  can 
assume  that  the  three  components  of  load  are  not  coupled.  A 
linearized  version  of  test  stage  position  calculation  is  also  used. 
These  algorithms  are  used  to  reduce  the  computation  time  so  that  a 
high  control  cycle  frequency  can  be  attained.  Control  cycle  frequency 
of  30  Hz  can  be  obtained  using  the  PDP  11/03  microcomputer.  With  the 
computer  interface,  variations  of  this  study  can  be  accomplished  with 
minimal  changes  in  the  software  of  the  system. 

3.4.2.  Load  Distribution 

The  findings  of  these  lumbar  creep  studies  suggest  that  an 
increasing  fraction  of  the  total  axial  load  applied  to  the 
intervertebral  joint  is  supported  by  the  intervertebral  disc  as  the 
joint  creeps.  This  increase  in  anterior  load  observed  at  0  degrees 
flexion-extension  angle,  with  zero  A/P  shear  is  probably  a  combined 
effect  from  sliding  of  the  facet  joints,  and  from  the  differences  in 
material  properties  between  the  intervertebral  disc  and  the  facet 
capsules.  The  results  suggest  that  the  facets  slide  along  their 
articular  surfaces,  passing  load  to  the  facet  capsules,  which  then 
relax. 

This  study  is  the  first  to  measure  the  time  dependent 
distribution  of  the  anterior-posterior  loads  in  the  spine.  Facet 
loads  of  25%  reported  by  Yang  and  King  [38]  and  20%  reported  by 
Nachemson  [21]  for  quasi-static  conditions  are  comparable  to  our 
result  of  23%  at  the  onset  of  loading.  Previously,  however,  there  was 
no  direct  evidence  of  how  axial  load  is  being  supported  by  the  facet 
joint.  Yang  and  King  [38],  Adams  and  Hutton  [1],  and  Andersson  [3] 
have  suggested  that  the  capsular  ligament  plays  an  important  role  in 
supporting  the  axial  compressive  load  across  the  facet  joint.  They 
suggested  that  the  ligament  is  visco-  elastic  and  relaxes  when  the 
deformation  of  the  capsular  ligament  is  held  constant.  As  each  facet 
slides  along  the  surface  of  the  adjacent  facet  in  response  to  an  axial 
load,  the  capsular  ligament  is  stretched,  so  as  to  resist  further 
sliding.  However,  the  relaxation  time  constant  of  the  capsular 
ligament  and  the  intervertebral  disc  are  different.  The  results  of 
our  experiments  suggest  that  the  relaxation  time  constant  of  the  facet 
is  shorter  than  that  of  the  intervertebral  disc.  This  indicates  that 
the  capsular  ligament  relaxes  faster  than  the  intervertebral  disc  when 
both  are  held  under  the  same  deformation.  The  long  and  short  time 
constant  calculated  from  the  Kelvin  solid  model  seem  to  agree  with 
this  hypothesis. 


3.4.3.  Relaxation  of  Lumbar  Specimens 


In  these  experiments,  the  spinal  specimens  were  tested  in  creep 
at  0  degrees  flexion-extension,  then  at  2  degrees  extension,  and 
finally  at  2  degrees  flexion.  Each  creep  test  lasted  for  1  hour  and 
the  specimens  were  allowed  to  relax  for  1  hour  between  tests.  Before 
relaxation,  the  specimens  were  positioned  at  the  original  neutral 
orientation  in  the  Planar  Testing  Apparatus. 

The  experimental  data  show  that  the  specimens  did  not  recover 
completely  from  the  creep  test  in  one  hour  and  that  the  disc  load 
level  of  the  2  degrees  extension  and  the  2  degrees  flexion  were  not 
consistent  between  tests,  Figure  24.  In  addition,  the  disc  load  at  2 
degrees  flexion  did  not  support  a  load  close  to  100%  of  the  total 
applied  axial  load  (1200  N)  as  was  suggested  by  Adams  and  Hutton  [1]. 
However,  the  first  test  performed  after  12  hours  of  relaxation  (0 
degrees  flexion-  extension)  always  produced  a  consistent  level  of  load 
starting  at  77%  of  the  total  applied  load  and  finished  at  89%  of  the 
total  applied  load,  Figure  25.  This  phenomena  was  confirmed  for  one 
specimen  when  the  displacement  history  was  repeated  on  the  intact 
specimen  (15].  The  results  show  that  the  load  history  (constant  at 
1200  N)  was  reproducible  after  the  intact  specimen  relaxed  for  12 
hours,  but  the  load  history  was  not  reproducible  if  the  specimen 
relaxed  for  only  2  hours. 

Koeller  has  reported  that  the  deformation  was  less  when  the 
specimen  was  creep  tested  a  second  time  [13].  However,  according  to 
Koeller,  the  difference  in  displacement  between  the  first  and  the 
second  test  was  very  small.  Figure  21  shows  that  the  displacement  of 
the  0  degrees  flexion-extension  test  is  consistently  the  lowest  and 
the  2  degrees  flexion  test  is  the  highest.  If  the  intervertebral 
specimens  are  not  given  enough  time  for  relaxation,  a  change  in 
viscoelastic  material  properties  of  the  intervertebral  joint  will 
occur.  Even  though  the  change  of  load  distribution  characteristic  due 
to  insufficient  relaxation  cannot  be  determined  from  these 
experiments,  the  data  shows  that  there  is  an  increase  in  deformation 
of  the  intervertebral  joint.  Based  on  these  findings,  the  relaxation 
time  between  tests  changes  the  load-displacement  characteristic  of  the 
spinal  specimens.  However,  more  data  are  needed  to  further  understand 
the  change  in  material  properties  of  the  intervertebral  joint  due  to 
repeated  creep  tests. 

3.4.4  Neutral  Orientation 

Earlier  stiffness  tests  of  the  intervertebral  joints  were  done 
using  constant  load  with  series  of  strings  and  pulleys  [17,25,26,31]. 
In  a  typical  stiffness  test,  the  specimen  is  displaced  by  a  constant 
static  load  and  measurement  of  the  displacement  made  after  5  to  10 
minutes.  Then  the  load  is  removed  and  the  specimen  is  allowed  to 
recover  for  10  to  15  minutes.  The  orientation  of  the  specimen  after 
the  recovery  period  is  defined  as  the  new  neutral  orientation.  During 
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relaxation,  the  specimen  may  or  may  not  return  to  its  original  height 
depending  on  the  specimen  conditions  such  as  disc  grade,  moisture,  and 
the  duration  of  the  experiment.  Hence,  the  new  neutral  orientation 
may  not  be  the  same  as  before. 

The  design  of  our  test  fixture  allows  creep  tests  to  start  at  the 
same  neutral  orientation  in  every  test  so  that  the  same  displacement 
condition  can  be  repeated  after  the  posterior  elements  are  transected. 
A  well  defined  neutral  orientation  is  important  because  a  shift  in  the 
neutral  orientation  can  lead  to  misintepretat ion  of  the  data 
especially  at  high  load  or  high  rates  of  loading  where  hysteresis  is 
pronounced  and  the  test  is  carried  out  beyond  the  toe  region  of  the 
load -d i splacemen t  curve. 


3.5.  SUMMARY  OF  LUMBAR  SPINE  STUDY 

These  experiments  suggest  several  important  features  of  the  load 
distribution  behavior  of  lumbar  intervertebral  joints. 

(1)  The  fraction  of  the  axial  load  supported  by  the  intervertebral 
disc  increases  with  time  when  a  constant  axial  compressive  load  is 
applied  to  the  intervertebral  joint.  At  0  degree  flexion-extension 
angle,  the  increase  in  the  intervertebral  disc  load  is  linear  with 
time.  The  initial  load  supported  by  the  disc  is  at  77%  of  the  total 
applied  load  and  increases  to  about  89%  in  one  hour.  At  2  degrees 
flexion,  the  intervertebral  disc  supports  nearly  all  of  the  total 
applied  axial  compressive  load.  At  2  degrees  extension,  the 
intervertebral  disc  load  increases  exponentially  with  time  and 
approaches  an  asymptotic  value.  The  final  load  supported  by  the  disc 
is  less  than  the  final  load  supported  by  the  disc  at  0  degree  flexion- 
extension  angle. 

(2)  Exactly  defined  neutral  orientation  and  sufficient  relaxation  time 
between  creep  tests  is  key  to  consistent  test  results.  For  a  1  hour 
creep  test,  12  hours  of  relaxation  at  4  degrees  Celsius  is  sufficient 
if  the  spinal  specimen  is  returned  to  its  original  geometric 
configuration  in  a  fixture.  Insufficient  relaxation  of  the  spinal 
specimen  causes  higher  than  normal  deformation  of  the  specimen.  The 
load-displacement  curve  of  an  insufficiently 

relaxed  specimen  demonstrates  a  greater  toe  region  and  greater  slope 
than  the  sufficiently  relaxed  specimen. 

(3)  Mechanisms  for  low  back  pain  cannot  be  directly  interpreted  from 
these  mechanical  tests.  However,  it  appears  that  at  2  degree 
extension,  the  tip  of  the  facet  joint  may  touch  the  lamina  of  the 
adjacent  vertebral  body. 

This  is  a  potential  mechanism  which  may  produce  pain  after  long 
periods  of  standing.  The  back  pain  elicited  after  prolong  erect 
sitting  may  be  caused  by  the  high  stress  exerted  on  the  intervertebral 
disc  when  load  is  transfered  to  the  disc  from  the  facet  joints. 
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The  material  presented  here  demonstrates  the  time  dependent  load 
distribution  characteristis  of  the  lumbar  intervertebral  joint.  The 
experiment  has  shown  that  time  is  an  important  factor  when  considering 
the  load  distribution  of  the  intervertebral  joint.  However,  more 
studies  of  this  aspect  of  spinal  biomechanics  are  needed  to  gain 
further  unders tanding  of  the  causes  of  low  back  pain. 


4.  VISCOELASTIC  AND  STRENGTH  PROPERTIES  OP  CERVICAL  SPINE 

The  study  of  spinal  mechanics  provides  clinically  useful 
information  for  the  analysis  of  movement  disorders  arising  from 
congenital  malformation  or  injury.  Detailed  data  defining  the 
kinematics,  the  relationship  between  loads  and  displacements,  the 
vertebral  strength,  and  the  modes  of  failure  of  the  normal  spine  are 
required  to  define  the  tolerance  of  the  cervical  spine  to  injury  and 
to  determine  when  and  if  surgery  or  other  treatment  is  necessary  to 
repair  an  abnormal  spine.  This  part  of  the  project  focused  on  the 
measurement  and  quantification  of  cervical  spinal  motion  under 
physiologic  loading  iji  vitro  in  both  non-destructive  and  failure 
testing  situations.  Testing  was  done  on  eight  cervical  spines  taken 
from  cadavers.  Cervical  spines  were  tested  nondestructi vely  to  gather 
load  vs.  displacement  data  and  destructively  in  a  combination  of 
flexion  and  compression. 


4.1  REVIEW  OF  THE  LITERATURE 

4.1.1  Anatomy 

The  human  spine  is  composed  of  the  cervical,  dorsal  (thoracic), 
lumbar,  sacral,  and  coccygeal  segments.  The  upper  segments  consist 
respectively  of  seven,  twelve,  and  five  separate  vertebrae,  while  the 
sacrum  and  the  coccyx  consist  respectively  of  five  and  four  vertebrae 
fused  in  the  adult.  The  spinal  column  has  four  curves  when  viewed 
laterally.  In  the  cervical  and  the  lumbar  segments,  the  spine  is 
convex  anteriorly,  and  in  the  thoracic  and  sacral  segments,  it  is 
concave  anteriorly.  Each  cervical,  thoracic,  and  lumbar  vertebrae 
consists  anteriorly  of  a  vertebral  body  and  posteriorly  of  a  neural 
arch  containing  four  articular,  two  transverse,  and  one  spinous 
process,  Figure  27.  The  arch  provides  a  protective  canal  for  the 
spinal  cord  and  nerves. 

Vertebrae  are  connected  to  one  another  anteriorly  through  the 
longitudinal  ligaments  and  the  f ibrocar t i lagenous  intervertebral 
discs,  and  posteriorly  through  the  interspinous  and  superspinous 
ligaments  and  the  facet  joints.  The  facet  joints  are  formed  by  the  two 
inferior  articular  processes  of  one  vertebra  and  the  two  superior 
articular  processes  of  the  vertebra  below  it.  The  transverse 
processes  project  laterally  and  serve  as  attachment  sites  for  the 
muscles  which  control  rotation  and  lateral  bending  of  the  spinal 


column.  The  spinous  process  projects  dorsally  and  distally  and  serves 
as  an  attachment  site  for  the  muscles  of  extension  and  rotation. 
Several  texts  present  a  thorough  description  of  the  anatomy  of  the 
normal  cervical  spine,  and  White  and  Panjabi  present  a  comprehensive 
description  of  the  biomechanics  of  the  entire  spine  [28). 

Spinal  movements  occur  between  motion  segments,  also  called 
functional  spinal  units  (FSU's),  consisting  of  two  adjacent  vertebrae 
and  their  interconnecting  tissue  and  ligaments.  The  FSU  is  the 
smallest  basic  spinal  structural  unit  which  exhibits  biomechanical 
characteristics  similar  to  the  entire  spine.  In  the  cervical  region 
the  seven  vertebrae  are  designated  Cl  through  C7.  The  cervical  spine 
consists  of  7  FSU's  beginning  with  articulation  between  the  head 
(occiput)  and  Cl,  Occ-Cl,  and  ending  with  C6-C7.  Motion  between 
vertebrae  takes  place  in  all  three  anatomic  planes,  the  sagittal, 
transverse  (horizontal),  and  coronal  (frontal)  planes.  The  four 
predominant  motions  of  spine  are  flexion-extension,  compression- 
tension,  axial  rotation,  and  lateral  bending. 

4.1.2  Biomechanical  Considerations 

Based  on  these  anatomic  considerations,  the  spine  is  an 
articulated  column  composed  of  structural  units  of  complex  geometries 
and  material  properties.  The  load  characteristics  (along  with  the 
multiaxial  nature  of  the  applied  loads)  result  in  additional 
complexities  at  the  structural  level.  These  complexities  of  spinal 
column  motion  are:  1)  coupled  motions,  2)  non-linearities  in  load 
versus  displacement  behavior,  and  3)  time  dependent  load-displacement 
responses.  These  were  recognized  in  previous  mechanical  studies  of 
connective  tissue  and  spinal  segments.  White  and  Panjabi  discuss  in 
detail  the  biomechanically  relevant  anatomy  [28]. 

Lysell  [19],  Panjabi  et  al  [25,26],  and  Markoff  [20]  and  others 
showed  that  spinal  motions  in  all  three  mobile  regions  are  generally 
coupled.  That  is,  motion  (translations  or  rotations)  in  one  anatomic 
direction  is  accompanied  by  motion  in  at  least  one  other  direction. 
For  example,  in  left  lateral  bending  of  the  distal  cervical  spine,  the 
vertebrae  also  rotate  axially  bringing  the  spinous  processes  away  from 
the  direction  of  bending.  Coupling  also  occurs  in  the  cervical  region 
between  flexion-extension  and  sagittal  anterior-posterior  translation 
[19].  The  coupling  of  spinal  motions  is  due  to  the  geometry  of  the 
vertebrae,  the  orientation  of  the  intervertebral  disc  and  facet 
articulation,  the  interconnecting  ligaments,  and  the  curvature  of  the 
spinal  column. 

The  load  vs.  displacement  response  of  the  spine  is  also  non¬ 
linear  [7,27,36].  When  forces  and  moments  are  applied  to  the  spine  and 
the  resulting  displacements  are  measured,  the  slope  of  the  load  vs. 
displacement  curve  increases  with  increasing  load,  Figure  28. 
Therefore,  displacement  measurements  made  in  uniaxial  loading 
situations  cannot  be  summed  in  order  to  indicate  the  response  during 
simultaneous  loading  along  more  than  one  axis  [7],  Combined  loading 
conditions  should  therefore  be  used  experimentally  whenever  possible. 


The  rate  of  loading  also  influences  the  load-deflection  behavior 
of  the  spine  [12,15,16].  This  is  a  result  of  the  viscoelastic 
characteristics  of  the  ligaments,  discs,  and  muscle.  The  rate  of 
displacement  is  directly  related  to  the  loads  produced  by  these 
components.  The  slope  of  the  load  vs.  displacement  curve  increases  at 
higher  loading  rates  and  the  forces  generated  are  generally  higher  at 
well.  Therefore,  a  thorough  study  of  spinal  motion  should  include 
several  loading  rates.  A  second  result  of  viscoelasticity  is  the 
hysteresis  evident  in  the  load-deflection  behavior  of  spines  loaded 
cyclically.  This  is  due  to  energy  loss  during  each  load  cycle,  Figure 
17. 


4.1.3  Objectives 

Detailed  information  on  the  load  vs.  displacement  response  of  the 
spine  can  best  be  gathered  through  in-vi tro  testing.  Past  spinal 
mechanics  investigators  have  generally  studied  two  vertebrae  segments 
(FSU's).  In  the  present  studies,  we  have  chosen  spinal  segments 
consisting  of  three  vertebrae  and  their  interconnecting  ligaments  and 
discs  (FVS's).  With  these  segments,  upper  and  lower  vertebrae  can  be 
manipulated  while  the  middle  vertebra  is  allowed  to  move  in  a  way 
physiologically  similar  to  the  iji  vivo  condition,  with  intervertebral 
discs  and  ligaments  being  the  only  constraints. 

Due  to  the  nonlinear  and  viscoelastic  nature  of  the  spine  and  the 
complex  loading  to  which  it  is  subjected,  a  meaningful  test  procedure 
must  include  a  multi  degree-of-f reedom  loading  apparatus  and  a 
displacement  measurement  device  which  does  not  interfere  with  the 
natural  modes  of  spinal  motion.  The  Planar  Testing  Apparatus, 
described  above,  is  capable  of  generating  either  controlled  loads  or 
controlled  displacements.  These  two  control  methods  make  possible  the 
measurement  of  all  of  the  parameters  necessary  to  define  independently 
specimen  stiffness  and  specimen  flexibility. 

The  displacement  measuring  apparatus  developed  for  this  project 
consists  of  three  precision  LVDT's  and  six  precision  RVDT's  grouped  so 
as  to  continuously  monitor  the  position  of  three  small  balls  on  a 
collar  rigidly  affixed  to  the  specimen.  The  collar  is  clamped  around 
the  free  vertebra  so  that  there  is  no  motion  of  the  collar  relative  to 
the  specimen.  However,  the  collar  is  designed  so  that  the  bone  is  not 
punctured,  thus  avoiding  unphysiologic  stresses.  The  transducer 
groups  are  designed  to  add  minimal  loading  to  the  specimen,  and  are 
capable  of  measurement  over  the  complete  range  of  normal  movement. 
The  displacement  collar  is  attached  only  to  the  moving  vertebra  and 
does  not  interfere  with  the  natural  motion  of  the  specimen.  The 
output  voltages  of  the  transducers  are  recorded  via  an  on-line 
computer.  With  this  apparatus,  a  highly  accurate,  continuous  and 
noninvasive  method  of  displacement  measuring  is  possible. 


The  motion  of  the  moving  vertebra  with  respect  to  the  fixed  can 
be  described  as  a  translation  along  and  Eulerian  rotation  about,  a 
fixed  orthogonal  axes  system  with  its  origin  at  the  center  of  the 
moving  body.  The  three  Euler  angles  very  nearly  describe  lateral 
bending,  flexion-extension,  and  torsion  rotations,  and  are  easily 
visualized,  particularly  for  planar  motion.  This  method  is  well 
suited  for  the  description  of  spine  kinematics. 

The  specific  objectives  of  this  cervical  spine  study  presented  in 
this  report  were: 

A)  Develop  experimental  techniques  needed  for  mechanical  testing  of 
human  cervical  spine  specimens. 

B)  Define  algorythyms  to  calculate  the  rigid  body  translations  and 
Euler  rotations  of  the  moving  vertebra  relative  to  a  fixed  reference. 

C)  Measure  the  nonlinear  mechanical  properties  of  intact  cervical 
spine  specimens  considering  sagittal  loading  and  displacements  for 
conditions  up  to  including  failure  of  the  vertebrae  and  rupture  of  the 
soft  interconnecting  tissues. 

A  final  objective  was  to  quantify  the  change  in  mechanical  properties 
and  strength  of  cervical  spine  segments  due  to  axial  rotation  applied 
in  combination  with  flexion  and  compression. 


4.2  MATERIALS  AND  METHODS 

4.2.1  Planar  Testing  Apparatus 

The  Planar  Testing  Apparatus  (PTA)  described  in  this  report  is  a 
multi  degree-of-f reedom  materials  testing  machine  capable  of 
multiaxial  load  or  displacement  control.  The  PTA  allows  in-vitro 

simulation  of  the  full  range  of  in-vivo  spinal  motions.  It  is  capable 
of  producing  simultaneous  multiaxial  displacements  and  loads  at 
constant  or  variable  rates.  System  software  presents  real  time  load 
vs.  displacement  data  for  monitoring  test  parameters.  While  in  load 
control,  direct  measurement  of  the  variables  required  to  determine 
spinal  segment  flexibility  can  be  made.  In  displacement  control, 
measurement  of  the  variables  necessary  to  determine  spinal  stiffness 
is  possible. 

4.2.2  Displacement  Measurement  Device 


A  complete  description  of  the  motion  of  a  single  vertebra  is 
defined  by  measurement  of  the  position  of  three  non-colinear  points 
fixed  to  the  moving  vertebra.  An  electromechanical  transducer  system 
was  developed  for  use  with  tie  PTA  for  the  measurement  of  the  absolute 
3-D  motion  of  the  center  vertebra  in  a  three  vertebrae  cervical  spine 
segment  ( 4 J .  The  3  vertebrae  segments  which  include  all  of  the 
interconnecting  ligaments  and  discs  will  be  referred  to  as  Functional 
Vertebral  Segments  (FVS's),  Figure  29. 
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This  six  degree  of  freedom  motion  measurement  system  consisted  of 
three  identical  transducer  groups  that  monitored  the  position  of  three 
points  attached  to  the  moving  vertebra.  Each  transducer  group  was 
composed  of  a  linear  variable  displacement  transducer  (LVDT)  and  two 
rotational  variable  displacement  transducers  (RVDT's)  (Schaevitz 
Engineering,  Cat.  No.  2000LV  and  R30A).  One  transducer  group  is 
illustrated  in  Figure  30.  The  tip  of  the  rod  from  the  transducer  was 
attached  to  a  point  on  a  collar  around  the  free  vertebra  being 
monitored.  The  transducer  rod  and  LVDT  housing  were  free  to  rotate 
about  the  axis  of  RVDT  1  and  RVDT  2.  RVDT  1  was  mounted  onto  the  LVDT 
housing  so  that  it  measured  motion  of  the  body  in  the  vertical  plane. 
RVDT  2  was  mounted  onto  the  housing  such  that  it  measured  motion  of 
the  LVDT  in  a  horizontal  plane.  The  LVDT  core  guide  was  attached  to 
the  rod  which  was  supported  in  two  precision,  linear  instrument 
bearings  housed  in  aluminum  blocks  attached  to  the  LVDT  body.  This 
arrangement  prevented  the  core  guide  from  being  stressed  during 
testing.  A  teflon  block  was  attached  to  the  back  of  the  LVDT  and  to 
the  supporting  rod  thus  allowing  free  translation  of  the  rod  along  the 
LVDT  body  with  no  lateral  displacement  or  rotation  of  the  core  of  the 
LVDT  in  the  LVDT  body.  The  output  from  the  RVDT's  was  linear  to 
within  0.25%  of  their  range  and  were  calibrated  with  a  sensitivity  of 
0.25  volts  per  degree  over  a  range  of  40  degrees.  The  LVDT's  were 
linear  to  within  0.25%  of  their  range  and  were  calibrated  with  a 
sensitivity  of  0.25  volts  per  mm  with  a  range  of  40  mm. 

Each  transducer  group  was  rigidly  attached  to  an  aluminum  bridge 
mounted  on  the  PTA  load  frame.  Figure  31  illustrates  the  orientation 
of  the  transducer  groups  on  the  bridge  with  respect  to  the  specimen 
and  the  load  frame.  The  groups  were  mechanically  balanced  so  that  no 
loads  were  applied  to  the  specimen  when  the  transducers  were  attached. 
With  the  nine  transducer  system,  it  was  possible  to  independently 
determine  the  3-D  location  of  three  non-colinear  points  attached  to 
the  specimen.  With  nine  transducers  monitoring  the  motion  of  a 
specimen  moving  with  6-D0F,  redundant  measurements  were  available  to 
check  accuracy  and  to  provide  back-up  measurements  in  case  of 
transducer  failure. 

4.2.3  Specimens 

Human  spines  were  collected  within  48  hours  of  death  from  Harvard 
Medical  School  through  its  Anatomical  Gifts  Program  (Table  7).  The 
spines  were  separated  into  lumbar,  thoracic,  and  cgrvical  segments, 
placed  in  plastic  bags,  sealed,  and  frozen  at  -20  C.  Freezing  and 
thawing  cycles  have  been  shown  by  others  not  to  affect  the  properties 
of  the  bone  [32],  the  disc  [11],  or  the  ligamentum  flavum  [2].  Befgre 
testing,  the  cervical  segments  were  allowed  to  thaw  overnight  at  5  C. 
Anterior-posterior  and  lateral  radiographs  were  then  taken.  Specimens 
showing  any  abnormalities,  fracture,  or  bone  disease  were  excluded. 
The  specimens  were  then  separated  by  dissection  at  the  discs  into 
three  functional  vertebral  segments  (FVS's),  either  C2-C3-C4  and  C5- 
C6-C7  segments  or  C3-C4-C5  and  C6-C7-T1  segments.  The  segments  were 
cleaned  of  excess  soft  tissue  while  leaving  the  ligaments  intact.  Two 
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steel  pins  1.75  mm  in  diameter  were  inserted  laterally  through  the 
inferior  and  superior  bodies  at  their  midlines.  The  pins  helped  to 
secure  the  bodies  in  a  castolytic  resin  ('Castolite'  Buehler  Ltd.  No. 
20-8120-009)  in  which  they  were  ultimately  embedded. 

4.2.4  Test  Fixture 

The  cervical  spine  holding  fixture  constructed  for  use  with  the 
PTA  consisted  of  an  aluminum  cup  114  mm  in  diameter,  12.7  mm  wall, 
rigidly  attached  to  a  25.4  mm  thick  aluminum  plate.  This  fixture 
assembly  was  rigidly  attached  to  the  test  stage  of  the  PTA.  Four  set¬ 
screws  9.5  mm  in  diameter  positioned  90  apart  were  screwed  through 
the  cup  wall  to  position  the  specimen  and  insure  that  the  resin 
remained  attached  to  the  cup.  Figure  32  illustrates  the  test  fixture. 
The  inferior  body  was  positioned  in  the  center  of  the  lower  cup  so 
that  the  plane  of  the  disc  was  parallel  to  the  cup  edge  and  to  the 
lower  plate.  Resin  was  then  poured  into  the  cup  covering  the  4  set 
screws,  2  pins,  and  approximately  3/4  of  the  height  of  the  body.  The 
2  exposed  vertebrae  were  wrapped  in  water  soaked  cotton  sheeting  and 
covered  with  saranite  film  to  prevent  dehydration  while  the  resin 
dried,  approximately  6  hours  at  5°  C.  This  technique  prevented 
specimen  dehydration  and  preserved  the  properties  of  the  disc  and  the 
longitudinal  ligaments  [11,37]. 

A  collar  made  of  galvanized  steel  was  used  for  attachment  of  the 
LVDT's  to  the  middle  vertebra,  Figure  33.  Three  3.9  mm  diameter  balls 
from  a  ball  and  cup  linkage  set  ('Du-Bro'  Cat.  N0.I8D)  were  screwed 
into  the  collar  in  a  radially  symmetric  pattern  to  provide  the  three 
reference  points.  The  cups  were  attached  to  the  tips  of  each  LVDT 
rod.  The  ball  linkages  provided  for  approximately  50  degrees  of 
rotation  with  very  little  friction.  Along  each  edge  of  the  collar 
were  cut  numerous  teeth-like  projections  about  3  mm  in  length  which 
held  securely  onto  the  vertebra  about  its  midline  without  puncturing 
the  bone.  The  collar  was  placed  around  the  body  and  fastened  in 
place.  Its  length  was  adjusted  to  fit  varying  sizes  of  specimens. 

Once  the  collar  with  its  3  balls  was  positioned,  the  specimen  was 
again  radiographed  laterally,  anterior-posteriorly ,  and  axially  so 
that  the  exact  3-D  orientation  of  the  balls  with  respect  to  the 
midpoint  of  the  free  vertebra  could  be  determined.  For  determining 
the  magnification  of  the  image,  2  small  pins  were  placed  into  the 
walls  of  the  inferior  aluminum  cup  at  about  the  midline.  The  actual 
distance  between  pins  was  compared  to  the  measured  distance  on  the 
image  to  determine  a  magnification  factor.  The  radiographs  were 
checked  to  see  that  the  collar  was  in  contact  only  with  the  center 
vertebra  of  the  specimen  and  its  posterior  elements.  The  superior 
body  was  then  positioned  and  embedded  in  a  second  cup-plate  fixture. 
The  upper  and  lower  plates  were  spaced  at  the  proper  distance  for 
embedding  the  specimen  and  held  parallel  through  the  use  of  aluminum 
tubes  of  the  required  height.  The  tubes  were  held  in  place  by 
threaded  rods  and  bolted  onto  the  upper  and  lower  plates.  Anterior- 
posterior  and  lateral  radiographs  were  taken  of  the  test  fixture  and 
the  specimen  so  that  the  orientation  of  the  center  of  the  body  of  the 


middle  vertebra  with  respect  to  the  test  fixture  could  be  determined. 
The  fixture  provided  a  method  for  securely  holding  the  specimen  so 
that  it  could  be  remounted  for  each  test  in  exactly  the  same  position. 
Figure  34  shows  a  collared  specimen  in  the  test  fixture  on  the  test 
stage  of  the  PTA. 

4.2.5  Test  Setup 

The  test  fixture  was  bolted  onto  the  test  stage  and  onto  the  load 
cell.  The  load  control  option  of  the  PTA  was  used  to  position  the 
test  stage  such  that  all  three  sagittal  loads  were  approximately  zero 
immediately  after  bolting  the  specimen  fixture  on  the  PTA.  The 
aluminum  tubes  and  the  threaded  rods  which  supported  the  upper  and 
lower  plates  of  the  fixture  while  mounting  the  specimens  were  then 
removed . 

Specimen  load-displacement  response  and  failure  test  data  were 
referenced  to  the  neutral  point  of  the  load-displacement  curves.  The 
neutral  point  is  defined  as  the  position  on  the  load-deflection  curve 
where  virtually  no  load  results  for  small  displacements.  The  neutral 
point  falls  within  the  flat  region  of  the  hysteresis  curve.  For  these 
tests,  the  neutral  point  was  defined  as  the  condition  midway  between 
the  two  zero  load  points  along  the  loading  curve  and  the  unloading 
curve  using  lmm/sec  and  1  deg/sec  displacement  input  rates.  At  this 
point,  there  were  no  forces  or  moments  applied  to  the  specimen.  Once 
the  neutral  point  had  been  determined  from  the  uniaxial  load- 
deflection  curves,  the  test  stage  was  moved  to  bring  the  specimen  to 
the  neutral  orientation.  Again,  the  test  stage  was  moved  in  small 
uniaxial  increments  and  the  resulting  load-deflection  curves  checked 
to  see  that  the  specimen  was  at  the  neutral  orientation.  The  test 
stage  location  was  then  recorded  so  that  the  specimen  could  be 
repositioned  at  the  original  neutral  orientation  for  future  testing. 
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geometric  center  of  the  moving  vertebral  body,  S  ,  was  referenced  to 
this  point.  The  loads  were  resolved  to,  and  the  moving  specimen 
rotated  about,  Sm.  This  point  is  at  the  intersection  of  the  midpoints 
of  ml  and  m2.  For  all  subsequent  testing,  the  PTA  was  able  to  resolve 
the  specimen  loads  to  the  neutral  point  and  to  generate  motions  of  the 
test  stage  with  reference  to  this  point. 

4.2.6  Test  Protocol 

Nondestructive  Testing 

Seven  specimens  were  tested  nondestructively  to  obtain  load  vs. 
displacement  behavior  at  two  different  displacement  rates.  Using  the 
load  control  option  of  the  PTA,  a  compressive  load  of  about  10£  body 
weight  was  first  applied  to  account  for  the  weight  of  the  head  and  the 
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stiffness  of  the  missing  muscles.  Hysteresis  curves  were  generated 
using  the  displacement  control  option  of  the  PTA.  The  specimen  was 
subjected  to  compression  and  tension  and  anterior  and  posterior  shear 
at  a  quasistatic  rate  of  1  mm  per  second  and  to  sagittal  bending  at  a 
rate  of  1  degree  per  second.  Each  test  was  repeated  three  times  and 
the  third  data  set  from  each  test  was  recorded.  This  was  to  ensure 
that  the  results  were  repeatable.  The  test  sequences  were  then 
repeated  at  a  rate  of  3  mm  per  second  and  3  degrees  per  second.  The 
third  data  set  from  each  test  was  recorded  to  determine  the  influence 
of  displacement  rate  on  the  load-displacement  response. 

Destructive  Testing 

After  completion  of  the  nondestructive  test  sequences,  each 
specimen  was  subjected  to  flexion  angles  of  up  to  42°  at  a  rate  of  6° 
per  second  until  failure  occured.  If  there  was  no  evidence  of 
ligamentous  or  bony  failure,  the  specimen  was  then  compressed  while 
held  at  maximum  flexion  until  failure  occurred.  The  forces  and 
moments  for  these  test  conditions  were  recorded  throughout  these 
tests.  Radiographs  were  taken  after  the  test  sequences  were  complete 
to  document  evidence  of  failure. 

4.2.7  Kinematic  Analysis 

A  FORTRAN  program  was  used  to  calculate  the  relative  motions 
between  the  fixed  and  the  moving  vertebra.  The  three  dimensional 
translation  vector  of  the  center  of  the  moving  vertebral  body,  S  , 
with  respect  to  the  center  of  the  fixed  vertebral  body,  S^,  and  tne 
Eulerian  rotations  of  the  moving  body  (Appendix  B)  were  calculated. 
The  translation  equations  and  the  rotation  equations  will  be  discussed 
separately . 

Translation 

The  3D  displacements  of  the  3  balls,  A,  B,  and  C,  on  the  collar 
were  determined  with  respect  to  the  pivot  points  of  their  respective 
transducer  groups,  A  ,  B  ,  and  C  .  The  transducer  voltages  were 
measured  and  transf8rmed°  into  distances  (mm)  and  angles  (radians) 
based  on  the  available  instrument  calibration  information  (Appendix 
A).  Figure  36  illustrates  the  vector  relationships  between  the  center 
of  the  fixed  and  moving  vertebra  (Sf,  S  ),  one  of  the  balls  attached 
to  the  moving  vertebra,  (A),  and  the  pivot  point  of  one  transducer 
group,  (Aq).  The  vector  (A-Aq)  can  be  determined  from 

(A-Aq)  =  [ZsJA  [Y]a  [Z]a  (L)a  (3) 

where 

(A-A  )  is  the  vector  in  the  pivot  reference  frame  from  the  point 
0  A  to  the  ball  A; 

[Zs)A  isa  coordinate  rotation  matrix  whose  columns  contains 

direction  cosines  of  the  transducer  group  axes  relative  to 
their  calibration  orientation; 
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[Y]  is  a  coordinate  rotation  matrix  whose  columns  contain  the 
direction  cosines  of  the  angle  determined  by  R1 ; 

[ z !  A  is  a  coordinate  rotation  matrix  whose  columns  contain  the 
rt  direction  cosines  of  the  angle  determined  by  R2 ;  and 

<l)a  is  the  vector  along  X,  to  the  point  A  whose  length  is 
determined  by  LI. 


The  elements  of  Equation  (3)  are  given  in  detail  in  Appendix  B. 

Similar  equations  relate  B  to  B  and  C  to  C  . 

oo 

The  positions  of  the  points  A  ,  B  ,  and  C  with  respect  to  the 
center  of  the  fixed  body,  S ^ ,  were  determined  trom  the  PTA  geometry 
and  from  radiographs  of  the  specimen  in  the  test  fixture.  The  vector 
position  of  ball  A  with  respect  to  was  then  determined  from 

(A)  =  <Aq)  *  (A-Aq)  (4) 

where 

(A)  is  the  vector  in  the  fixed  reference  frame,  Sp  from  the 
origin  of  the  fixed  frame  to  the  point  A;  and1 

(Aq)  is  the  vector  in  the  fixed  reference  frame,  Sp  from  the 
origin  of  the  fixed  frame  to  the  point  Aq. 

The  position  vectors  from  ball  A  to  and  from  ball  A  to  were 
determined  from  radiographs.  The  equation  for  the  position  vector 
from  S  to  Sf  is 

m  (S)  =  (A)  H-  (R](S)a  (5) 

where 

(S)  is  the  position  vector  from  the  origin  of  the  moving 
reference  frame,  (S  )  to  the  origin  of  the  fixed 
frame,  (S^);  m 

(A)  is  the  position  vector  in  the  fixed  reference  frame  from  the 
origin  of  the  fixed  frame,  Sp  to  the  ball  A; 

[RJ  is  the  rotation  matrix  whose  elements  represent  successive 

rotations  of  S  about  X  ,Y  ,  and  Z  ;  and 
m  m  m  m 

(S^)  is  the  position  vector,  in  the  balls  reference  frame,  of  Sm 
with  respect  to  A.  The  elements  of  were  determined  from 
radiographs. 


Eulerian  Rotations 

Figures  36  and  37  illustrate  the  vectors  used  to  calculate  the 
unknown  Euler  rotations  of  the  point  S  .  A  method  was  chosen  that 
simultaneously  determines  both  the  three  translations  of  S  with 
respect  to  the  fixed  coordinate  system,  Xp  Y^,  Z , ,  and  the  successive 
Eulerian  rotations  about  the  fixed  axes.  This  method  was  used 
previously  by  Dieudonne,  et  al.  at  NASA  [6]  to  determine  the  position 
and  orientation  of  a  payload  platform  of  a  motion  simulator  relative 
to  a  fixed  platform,  and  used  by  Koogle,  et  al.  [14]  to  describe  the 
motions  of  a  special  transducer  system.  The  resulting  translations 
were  compared  with  those  determined  in  the  previous  section  to  check 
on  the  accuracy  of  the  methods.  The  nomenclature  presented  by 
Dieudonne,  et  al.  was  used  here. 


X  ,  Y  ,  Z  mutually  perpendicular  set  of  unit  vectors  defining 
the  moving  reference  frame  with  its  origin  at  the 
center  of  the  moving  vertebral  body,  S  ; 

X  j ,  'i  r ,  mutually  perpendicular  set  of  unit  vectors  defining 
1  the  fixed  reference  frame  with  its  origin  at  the 

center  of  the  fixed  vertebral  body,  S^; 


A.  vector  fixed  in  the  moving  reference  frame  from  S  to  the 
1  balls  A,  B,  and  C;i  =  l,2...6;  11 

A.  ,  A.  ,  A.  components  of  vector  A.; 

ix  ly  iz  l 

B.  vector  fixed  in  the  fixed  reference  frame  from  S,  to  points 
1  A' ,  B' ,  and  C'  on  the  PTA  frame;  i  =  1 , 2 . . . 6 ; 

B.  ,  B.  ,  B.  components  of  the  vector  B.; 
ix’  iy’  iz  K  i 

S  vector  in  the  fixed  reference  frame  from  S£  to  S  ; 

t  m 

x,  y,  z  components  of  the  vector  S; 

1.  vector  in  the  fixed  reference  frame  from  points  A',  B' ,  and  C' 
to  the  points  A,  B,  and  C;  i = 1 , 2 . . . 6 ; 

1.  ,  1.  ,  1.  components  of  the  vector  1.; 

|1*|  magnitude  of  1.  representing  the  actual  length  as 
1  3  determined  from  the  PTA  geometry  and  the  locations  of  A,  B, 
and  C; 

\p,  9,  component  rotations  of  the  moving  body  about  the  X^,  Y^ 
and  Zj  axes  respectively; 

[R]  rotation  matrix  representing  component  rotations  of  the 
moving  body  relative  to  the  fixed  body; 

f(a)  vector  function; 
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0  null  vector; 
a  unknown  parameter  vector. 


From  the  transducer  outputs,  the  components  of  the  vectors  1.  are 
known.  This  vector  can  also  be  given  by  1 


1.  =  [Rj  A.  +  S  -  B.,  i  =  1 , 2 . . .6.  (6) 

Determination  of  the  relative  orientation  of  S  with  respect  to  Sf 
involves  solving  six  simultaneous  nonlinear  equations  represented  by 
Equation  6  for  the  six  unknowns  (x,  y,  z,  <p,  9,  $) .  The  Newton-Raphson 
iteration  method  was  used  to  compute  the  vector  root,  a,  of  the 
equation 


f(a)  =  0. 


This  function  is  defined  as 


f .  ( a)  =  l.Tl.  -  |1.|2  . 

The  iteration  formula  is  1  1  1  a 


a  ,  =  a  - 
n+1  n 


8f(an) 


-1 


£(-n> 


where 


and 


(7) 


(8) 

(9) 


(10) 


3f(an) 

-  and  f(a  )  defined  in  detail  in  Appendix  B. 

8a  n 

n 


with 
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Given  the  six  1.  values  and  initial  guesses  for  x,  y,  z,  0,  and  4>, 
Equation  6  is  used  iteratively  until  a  given  convergence  criteria  is 
reached  which  defines  the  position  and  orientation  of  the  center  of 
the  moving  body,  S  relative  to  the  center  of  the  fixed  body,  S^.  The 
convergence  criterion  was 

t  r  I  a  ,  -  a  I  <  0 . 0001 

1  n+1  n  1  - 

where  tr(a)  is  the  sum  of  the  components  of  a. 


4.3  RESULTS 

The  results  are  presented  in  the  form  of  load  vs.  displacement 
curves  made  from  a)  the  data  obtained  from  the  load  cell,  b)  the 
transducers  monitoring  the  motions  of  the  test  stage,  and  c)  the 
specimen  collar.  The  curves  for  all  specimens  are  plotted  using  the 
same  scale.  An  expanded  form  of  each  curve  is  presented  in  Appendix 
D.  The  three  components  of  load  (anterior-posterior  shear,  tension- 
compression  force,  and  flexion-extension  moment)  were  resolved  from 
the  load  cell  to  the  center  of  the  middle  vertebral  body,  and  are 
plotted  against  the  displacements  of  the  middle  vertebra  with  respect 
to  the  center  of  the  fixed  vertebra  above  it.  Data  were  obtained 
throughout  the  nondestructive  test  sequences  and  for  the  destructive 
test  sequences  when  angles  of  rotation  were  below  15°.  For  higher 
angles  of  rotation,  the  displacement  measuring  collar  was  removed. 
For  destructive  test  sequences,  the  loads  are  plotted  against  test 
Stage  displacement. 

The  data  analysis  focuses  on  two  biomechanical  considerations: 
coupling  of  spinal  motions  and  the  effect  of  displacement  rates  on 
stiffness.  The  data  are  divided  into  two  groups,  one  obtained  from 
the  upper  cervical  region,  (specimens  C,  F,  H,  and  E),  and  the  other 
for  lower  cervical  region,  (specimens  B,  D,  and  G). 


4.3.1  Nondestructive  Axial  Tests 

The  results  of  the  compressive  force  vs.  axial  displacements  are 
shown  for  the  two  levels  of  the  cervical  spine  in  Figures  38.1-38.4 
(for  the  upper  cervical  region)  and  in  Figures  38.5-38.7  (for  the 
lower  cervical  region).  The  three  vertebrae  specimens  were  loaded  at 
compressive  rates  between  0.5  mm/sec  to  5  mm/sec.  For  those  specimens 
loaded  at  more  than  one  rate,  the  results  are  presented  on  the  same 
curve . 

The  compressive  response  of  all  specimens  displays  a  large  region 
of  low  load  with  negligible  slope  near  the  neutral  point.  The  slope 
increases  beyond  the  neutral  zone  in  both  tension  and  compression. 
Compressive  forces  are  shown  positive.  The  slopes  of  the  curves  from 
0.5  mm  to  0.7  mm  in  tension  and  0.4  to  0.6  mm  in  compression  were  used 
to  estimate  the  specimen  axial  stiffness.  If  more  than  one 


displacement  input  rate  was  tested,  the  results  from  the  lower  rate 
curve  were  used.  From  these  measurements,  all  specimens  were  stiffer 
in  compression  than  in  tension.  In  the  specimens  taken  from  the  upper 
cervical  region,  the  estimated  stiffnesses  ranged  from  145  to  325  N/mm 
for  tension,  and  from  500  to  3500  N/mm  for  compression  (Table  8).  Two 
of  these  compression  specimens  (F  and  H)  exhibited  very  high  stiffness 
values.  No  specific  reason  for  this  behavior  could  be  identified,  but 
bone  on  bone  contact  could  have  been  a  possibility.  In  the  specimens 
taken  from  the  lower  level  of  the  spine,  the  estimated  stiffness 
ranged  from  25  to  35  N/mm  in  tension  and  from  200  to  360  N/mm  in 
compression.  There  was  no  difference  in  the  slope  or  position  of 
these  curves  as  the  rate  of  loading  was  increased. 

Figures  39.1-7  compare  coupled  response  of  the  vertical  and 
shear  forces  vs.  vertical  displacement.  The  curves  of  shear  force  vs. 
Z-displacement  display  a  large  region  of  low  load  and  zero  slope  near 
the  neutral  point.  When  the  stage  was  displaced  along  Z  (compression- 
tension),  the  shear  forces  were  zero  in  tension  but  there  was  slight 
posterior  shear  associated  with  compression. 

Figures  40.1-7  compare  displacements  of  the  center  vertebra  along 
X,  Y,  and  Z  with  the  inputs  of  the  test  stage  displacement  along  Z. 
There  were  only  slight  displacements  along  Y  in  tension  and 
compression  for  all  specimens,  however  there  was  relatively  large 
posterior  displacement  associated  with  tension  for  specimens  C,  F,  and 
G  and  anterior  displacement  associated  with  compression  for  specimens 
C  and  F. 


4.3.2  Nondestructive  Shear  Tests 

Figures  41.1-7  show  the  measured  shear  force  plotted  against  the 
anterior-posterior  X-displacements.  These  force  vs.  displacement 
trends  are  similar  to  those  in  the  axial  force  vs.  Z-displacement 
plots.  The  slopes  were  near  zero  for  small  displacements,  and 
increased  as  anterior  or  posterior  displacement  increased.  Estimates 
of  the  stiffness  of  the  specimens  was  measured  from  the  slope  of  the 
curves  between  0.5  mm  to  0.7  mm  of  anterior  motion  and  0.6  mm  to  0.8 
mm  of  posterior  motion.  For  the  specimens  taken  from  the  lower  level 
of  the  spine,  the  anterior  stiffness  (25  N/mm-30  N/mm)  was  greater 
than  the  posterior  stiffness  (5  N/mm-14  N/mm),  although  all 
stiffnesses  were  fairly  small.  In  the  specimens  taken  from  the  upper 
levels,  specimens  F,  H,  and  E  dis  layed  much  higher  posterior 
stiffnesses  than  anterior  stiffnesses  The  estimated  stiffnesses 
ranged  from  30  N/mm  to  70  N/mm  for  anterior  shear,  and  from  35  N/mm  to 
250  N/mm  for  posterior  shear  in  the  specimens  from  the  upper  levels. 
The  slope  of  these  curves  did  not  suggest  a  dependence  of  estimated 
stiffness  on  displacement  rate  for  the  rates  tested. 

The  axial  and  shear  forces  generated  by  X-displacements  are 
compared  in  Figures  42.1-7.  Vhen  the  stage  was  displaced  posteriorly 
along  X,  the  compressive  forces  generated  were  approximately  50X 
greater  than  those  generated  in  shear  except  for  specimens  B  and  D. 
When  the  stage  was  displaced  anteriorly,  there  were  little  or  no 
vertical  forces  generated. 
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The  curves  in  Figures  43.1-7  show  the  resulting  rotations  for 
anterior-posterior  displacement  inputs.  When  the  lover  vertebra  and 
test  stage  was  displaced  posteriorly  along  X,  the  middle  vertebra  of 
each  specimen  rotated  positively  about  Y,  Figure  29.  When  the  stage 
was  displaced  anteriorly  along  X,  the  vertebra  rotated  negatively 
about  Y.  There  was  no  region  of  low  moment  and  zero  slope  near  the 
neutral  point  of  the  curves.  Specimens  G,  E,  and  H  show  some  lateral 
bending  associated  with  anterior-posterior  shear. 

4.3.3  Nondestructive  Flexion-Extension  Tests 

Figures  44.1-7  show  the  flexion-extension  moment  plotted  against 
sagittal  rotation  curves  for  each  specimen.  The  flexion  moment  and 
rotation  are  shown  negative.  The  results  of  the  bending  tests  show 
significant  scatter  at  the  low  angles  of  rotation  shown.  For  moment 
vs.  rotation  data,  the  stage  was  rotated  1.5°  to  6°  in  flexion- 
extension.  Based  on  the  data  collected  by  White  and  Panjabi  [ 38 ]  this 
is  about  15%  of  the  range  of  motion  of  the  typical  cervical  FVS.  In 
general,  the  specimens  were  stiffer  in  extension  than  in  flexion. 

The  axial  and  shear  forces  generated  when  the  specimen  was 
rotated  in  flexion-extension  are  shown  in  the  curves  in  Figures  45.1- 
7.  A  small  amount  of  compression  resulted  with  extension  rotation  and 
tension  was  measured  with  flexion,  except  for  specimen  H  which  was 
under  compression  in  both  flexion  and  extension. 

4.3.4  Destructive  Tests 

Figures  46.1-6  show  the  sagittal  moment  resulting  from  large 
flexion  angle  inputs.  A  modified  test  protocol  was  used  for  two  of 
the  specimens,  C  and  E.  Specimen  C  was  not  tested  at  large  angles  of 
flexion.  Specimen  E  failed  after  the  initial  nondestructive  testing 
was  complete  and  was  then  tested  at  large  angles  of  flexion  for 
comparison  purposes.  The  load-displacement  data  for  specimen  E  at 
the  time  of  failure  are  not  available.  Specimen  E  displayed  only  -800 
N-mm  of  moment  for  a  sagittal  rotation  of  -30°. 

Ag  shown,  for  example  in  Figure  46.1,  the  neutral  zone  extends 
from  0°  to  at  least  10°  before  the  slope  of  the  curve  increases.  When 
a  rotation  did  not  produce  evidence  of  failure,  the  specimen  was 
returned  to  its  neutral  position.  The  stage  was  then  rotated  to  a 
higher  angle.  The  curves  for  each  test  are  shown  on  the  figure.  This 
procedure  was  followed  for  specimens  F,  G,  and  H,  and  Figures  47.1-3. 
The  curves  showing  the  moment  resulting  from  each  flexion  sequence  are 
also  shown.  Each  time  the  test  was  repeated  to  a  higher  angle,  the 
slope  of  the  load-displacement  curve  decreased.  There  was  no  clear 
evidence  of  failure  on  the  moment  vs.  sagittal  rotation  curves  for  any 
specimen . 

The  curves  in  Figure  48.1-6  show  the  vertical  and  shear  forces 
generated  for  large  sagittal  rotations.  There  was  a  slight  tensile 
force  generated  for  low  angles  of  flexion.  This  axial  force  became 
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compressive  and  increased  markedly  as  the  angle  of  flexion  was 
increased,  however  the  force  remained  far  below  the  maximum 
compressive  force  of  1750  N  presented  by  Panjabi,  et  al.  [28]  for 
cervical  spine  from  C3-C7.  A  negative  shear  force  was  generated  as 
the  angle  of  flexion  increases.  Specimen  E,  which  was  broken  before 
these  tests,  showed  the  opposite  trends. 

If  there  was  no  evidence  of  yield  on  the  moment  vs.  sagittal 
rotation  curves  in  the  destructive  tests,  the  stage  was  rotated  to  the 
maximum  flexion  position  possible  with  the  test  fixture  and  then 
compressed.  Yield  on  the  load-displacement  curves  was  measured  at  the 
point  of  zero  slope,  maximum  load,  however  only  specimen  D  displayed 
clear  evidence  of  failure  on  the  curve  showing  vertical  force  vs. 
vertical  displacements,  Figure  49.  This  occured  at  an  angle  of  -27 
deg  and  6.5  mm  of  compression,  resulting  in  a  compressive  force  of 
2600  N. 

All  specimens  were  examined  visually  to  determine  which  elements 
had  failed  when  he  test  sequences  were  complete.  In  all  cases,  there 
was  complete  failure  of  posterior  ligaments  and  both  discs  adjacent  to 
the  middle  vertebral  body.  The  specimens  were  then  radiographed  to 
determine  if  bony  failure  had  occured.  The  radiographs  of  specimens  D 
and  H  showed  anterior  fracture  of  the  middle  vertebral  body  as  well  as 
complete  posterior  element  failure.  The  radiograph  of  specimen  G 
showed  fracture  of  the  anterior  superior  vertebral  body,  fracture  of 
the  superior  part  of  the  middle  vertebral  body,  and  complete  posterior 
element  fracture.  The  radiographs  for  the  remaining  three  specimens 
showed  no  bony  fracture. 

4.4  DISCUSSION 

4.4.1  Test  Results-Cervical  Segments 

All  specimens  show  greater  stiffness  in  compression  than  in 
tension  based  on  estimates  of  the  slopes  of  the  load  displacement 
curves.  This  was  typical  of  segments  from  both  regions  of  the  spine. 
This  is  apparently  due  to  the  compression  resistance  offered  by  the 
discs  and  the  compressive  stiffness  of  the  facet  joints.  Specimens  F 
and  H  showed  extremely  high  compressive  stiffnesses.  At  only  0.6  mm 
of  vertical  compression,  these  specimens  were  subjected  to  the  maximum 
force  that  cervical  vertebrae  can  withstand,  according  to  White  and 
Panjabi  [37],  yet  they  showed  no  obvious  signs  of  failure.  The 
reasons  for  this  are  not  apparent,  however  bony  interaction  between 
the  adjacent  vertebra  in  the  most  plausible  explanation. 

The  vertical  (compression-tension)  stiffnesses  decrease  from  the 
upper  cervical  levels  to  the  lower  cervical  levels.  This  trend  is 
evident  even  without  considering  specimens  F  and  H.  Liu,  et  al. 
[18]  studied  the  biomechanics  of  cervical  FSU's  at  rates  of  0.0085 
mm/sec  and  85  mm/sec.  They  compressed  a  FSU  steadily  until  a  load  of 
880  N  was  reached  and  then  returned  the  specimen  to  its  original 
configuration  at  the  same  rate.  In  tension,  a  maximum  displacement  of 
1.25  mm  was  recorded.  Liu,  et  al.  reported  that  axial  stiffness 
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increases  in  more  distal  segments  of  the  spine.  To  determine  how  his 
stiffness  data  from  FSU's  compared  stiffnesses  of  FVS's,  we  summed  the 
stiffnesses  from  two  adjacent  FSU's  as  if  they  were  springs  in  series 
(k^k^/k^+l^) •  The  results  of  this  calculation  showed  no  evidence  of 
increasing  or  decreasing  stiffness  in  the  cephalocaudal  direction. 

The  displacement  rates  used  in  the  nondestructive  test  sequences 
here  were  relatively  low  and  can  all  be  considered  quasistatic.  This 
explains  why  there  was  no  dependence  of  estimated  stiffness  on 
displacement  rate  for  either  axial  or  anterior-posterior  loading. 
Tests  conducted  at  rates  from  0.5  mm/sec  to  5  mm/sec  showed  no 
significant  changes  in  the  load-displacement  response.  Dynamic 
testing  is  generally  done  at  rates  approaching  100  mm/sec  to  simulate 
pilot  ejection  or  automobile  crash  environments  and  generally  show  an 
increase  in  the  slope  of  load-displacement  curves.  While  the  PTA  is 
capable  of  rapid  displacements,  we  were  not  able  to  take  high  speed 
data  during  these  tests  due  to  the  data  limitations  of  the  MINC 
computer . 

Shear  stiffnesses  were  consistently  lower  than  axial  stiffnesses. 
This  finding  is  typical  also  of  the  thoracic  and  lumbar  regions.  The 
compressive  forces  are  generated  by  the  interaction  of  the  facet 
joints  and  are  generally  larger  than  the  shear  forces  generated  by 
negative  X-displacements .  As  the  middle  vertebra  is  displaced 
anteriorly  with  respect  to  the  fixed  vertebra,  little  or  no  vertical 
forces  are  generated. 

In  the  specimens  from  the  upper  cervical  levels,  stiffnesses  were 
higher  in  posterior  shear  than  in  anterior  shear,  with  the  exception 
of  specimen  C.  This  is  due  to  the  interaction  of  the  apophyseal  joint 
members  as  the  middle  vertebra  is  moved  posteriorly  with  respect  to 
the  fixed  vertebra  above  it.  This  trend  is  reversed  in  the  specimens 
from  the  lower  levels,  although  the  stiffnesses  displayed  in  both 
anterior  and  posterior  shear  in  the  lower  segments  were  very  low. 
Liu,  et  al.  [18]  reported  that  posterior  shear  stiffnesses  were  higher 
than  anterior  shear  stiffnesses  for  all  levels  of  the  cervical  spine. 

The  shear  stiffness  also  decreases  from  the  upper  cervical  levels 
to  the  lower  levels,  with  the  exception  of  specimen  C.  This  may  be 
due  to  the  changing  geometry  of  the  more  distal  posterior  elements. 
Liu,  et  al.  [18]  reported  the  opposite  trend.  When  his  data  from  two 
adjacent  FSU's  are  summed  (k. k2/k^+k2) ,  there  does  not  seem  to  be  any 
evidence  of  increasing  or  aecreasing  stiffness  as  one  descends  the 
cervical  spine. 

The  middle  vertebra  rotated  positively  about  Y  (extension)  as  the 
vertebra  below  was  displaced  in  the  negative  X-direction,  again  due  to 
the  interaction  of  the  facet  joints  and  to  the  stiffness  of  the 
anterior  ligament.  It  rotated  negatively  about  Y  as  the  vertebra 
below  was  displaced  in  the  +X-direction  due  to  the  anterior 
compression  of  the  discs  and  to  resistance  offered  by  the  posterior 
ligaments.  The  coupling  of  horizontal  displacements  and  sagittal 
rotations  of  cervical  spinal  segments  has  been  shown  by  Liu,  et  al. 


[18]  and  Panjabi,  et  al.  [27].  Panjabi  applied  a  horizontal  load  of 
25%  body  weight  to  cervical  FSU's  and  measured  the  resulting  sagittal 
rotation.  He  reported  maximum  X-displacements  of  2.7  mm  and  maximum 
Y-rotations  of  10.7°.  In  the  present  study,  a  compressive  preload  of 
10%  body  weight  was  added  to  the  FVS  and  a  maximum  of  5U  extension 
resulted  from  a  maximum  of  1.5  mm  of  displacement  of  the  lower 
vertebra  posteriorly,  in  the  negative  X-direction  (Fig.  20.3).  The 
amount  of  preload  may  greatly  affect  the  degree  of  rotation  associated 
with  horizontal  (X)  displacements.  There  was  a  small  amount  of  axial 
rotation  about  Z  also,  most  notably  for  specimens  E,  G,  and  H.  This 
was  due  most  probably  to  irregularities  in  these  specimens.  Specimens 
G  and  H  are  from  the  same  cervical  spine. 

The  widest  scatter  in  our  data  between  specimens  was  measured  in 
the  curves  of  sagittal  moment  vs.  sagittal  rotation  (Fig.  21.1-21.7). 
The  angles  of  rotation  are  too  low  to  detect  any  trends.  In 
quasistatic  A-P  bending  tests  on  cervical  FSU's,  Liu,  et  al.  [18]  also 
found  bending  moments  to  vary  widely.  For  flexion  and  extension 
angles  of  5°,  he  found  maximum  extension  moments  ranging  from  400  N-mm 
to  7500  N-mm  and  maximum  flexion  moments  ranging  from  -800  N-mm  to  - 
9000  N-mm  with  no  dependence  on  cervical  spinal  level.  In  future 
testing,  the  angles  of  rotation  should  be  increased  to  at  least  30%  of 
the  physiologic  range  of  motion  of  cervical  FVS's  to  provide  better 
description  of  the  flexion-extension  properties. 

There  were  compressive  forces  generated  at  the  middle  vertebra 
when  the  specimens  were  rotated  in  extension  due  to  the  interplay  of 
the  facet  joints  and  the  intervertebral  discs.  Compressive  forces 
were  less  than  80  N  for  angles  of  extension  of  less  than  1.5°,  Figure 
45.2. 

There  is  no  evidence  of  sudden  failure  of  FVS's  at  high  angles  of 
flexion.  This  is  clear  from  the  curves  showing  moment  vs.  stage 
rotation,  Figures  46.1-6.  When  comparing  the  forces  generated  in  X 
and  in  Z  for  these  large  flexion  angles,  Figures  48.1-6,  it  is  clear 
that  these  forces  alone  cannot  account  for  failure  of  the  specimens. 
However,  the  study  of  Figures  47.1-3  suggest  that  each  time  the 
specimens  are  rotated  to  successively  higher  flexion  angles,  the  soft 
tissues  and  the  posterior  ligaments  become  less  stiff.  In  this  way, 
the  failure  of  the  specimen  may  be  a  gradual  giving  way  of  many 
structural  elements  at  once. 

When  specimen  D  was  compressed  while  held  at  27°  of  flexion, 
there  was  a  bony  fracture  at  2600  N.  Liu,  et  al.  [18]  found  a  maximum 
compressive  strength  (without  flexion)  of  1020  N  for  a  FSU  consisting 
of  C6-C7. 


4.4.2  Displacement  Measurement 

The  PTA  allows  simulation  of  the  full  physiologic  range  of  spinal 
motions.  This  machine  provides  continuous  and  repeatable  loading  or 
displacement  capabilities.  Together  with  the  displacement  measuring 
collar,  this  makes  possible  the  measurement  of  all  necessary 
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parameters  for  the  definition  of  the  complete  biomechanical  response 
of  the  spine,  including  both  specimen  stiffness  and  flexibility  and 
low  rate  viscoelastic  properties.  The  displacement  collar  provided  an 
ideal  means  of  monitoring  motion.  It  held  firmly  yet  non-invasively 
around  the  middle  vertebra. 

Electromechanical  transducer  assemblies  have  the  advantage  over 
dial  gage  systems  of  highly  accurate  and  continuous  real-time 
recording  of  motion.  Previous  investigators  have  attached 
electromechanical  devices  or  balls  to  the  moving  member  of  a  FSU  with 
bone  screws  or  Steinmann  pins.  While  these  methods  are  valid  for  non¬ 
destructive  testing,  they  would  most  likely  affect  the  outcome  of 
failure  testing  due  to  increased  stresses  at  the  site  of  the  screw  or 
pin.  Further,  Panjabi,  et  al.  [28]  noted  difficulty  in  keeping  his 
transducers  in  contact  with  the  balls  attached  to  the  moving  member. 
He  overcame  this  difficulty  by  adding  rubber  bands,  or  by  adding  low 
pressure  compressed  air  to  the  back  of  the  LVDT  coils.  This,  together 
with  the  weight  of  the  relatively  large  balls,  could  apply  additional 
stresses  to  the  specimen. 

The  six  degree  of  freedom  displacement  measuring  apparatus  used 
here  combines  the  accuracy  of  precision  transducers  with  a  lightweight 
and  virtually  friction  free  gimbaled  assemblies.  This  system  provides 
a  non-invasive  yet  rigid  attachment  to  the  specimen  needed  for 
accurate  measurements  without  hampering  normal  motions.  Since  a  total 
of  nine  transducers  are  used,  three  more  transducers  than  are 
necessary  to  determine  the  spatial  location  of  three  noncolinear 
points,  there  is  a  backup  system  in  the  event  that  one  or  more 
transducer  becomes  non-operat ional  during  testing. 

A  trial  displacement  collar  was  used  to  determine  the  accuracy  of 
the  entire  displacement  measurement  system  involved  the  use  of.  By 
using  the  collar  in  conjunction  with  the  PTA,  it  was  possible  to 
determine  the  errors  from  this  apparatus  for  a  situation  identical  to 
specimen  tests.  Curves  were  plotted  which  show  displacement  collar 
motion  and  test  stage  motion  vs.  time  (Appendix  C).  Minimum  errors  of 
close  to  OX  were  recorded  for  X  and  Y  translation  directions.  Minimum 
errors  of  3%  were  recorded  for  measurement  of  rotation  about  £he  X  and 
Z  axes.  Maximum  errors  of  1.2°  for  a  stage  rotation  of  10  about  Y 
and  0.4  mm  for  Z  stage  translations  of  3  mm  were  recorded.  In 
addition  to  transducer  error,  other  possible  sources  for  this  error 
included  motion  of  the  displacement  collar,  position  of  the  assembly 
on  the  PTA  bridge,  or  in  the  determination  of  the  position  of  the 
assembly  origins  with  respect  to  the  global  reference  system. 

4.5  SUMMARY  OF  CERVICAL  SPINE  STUDY 

This  study  of  the  cervical  spine  measured  the  biomechanical  load- 
displacement  properties  of  specimens  up  to  failure  using  a  combination 
of  refined  experimental  and  analytical  methods.  Previous 
biomechanical  investigations  of  spine  segments  and  the  current 
investigation  have  shown  that  the  spine  displays  complex  nonlinear, 
coupled,  viscoelastic  properties.  Therefore,  to  best  examine  this 


response,  physiological  load-displacement  conditions  must  be 
reproduced  in  vitro,  and  the  resulting  three  dimensional  motions 
measured  and  analyzed. 

The  use  of  three-vertebrae  Functional  Vertebral  Segments  (FVS)  in 
mechanical  tests  offered  advantages  over  both  two  vertebrae  segments 
(FSU's)  and  whole  cervical  spines.  Tests  of  the  FSU  involve  fixing 
one  vertebra  and  applying  loads  to  the  second,  taking  into  account  the 
motions  about  one  joint  only.  Accurate  strength  measurements  cannot 
be  obtained  from  FSU's.  For  the  FVS,  the  middle  vertebra  is  jointed 
naturally  to  the  two  adjacent  vertebrae,  providing  for  more 
physiological  end  conditions.  Unlike  whole  cervical  spine  specimens, 
motions  of  each  vertebra  can  be  controlled  and  directly  measured.  The 
use  of  the  Functional  Vertebral  Segment  is  recommended  in  future 
studies  which  include  failure  investigations. 

A  valid  test  apparatus  for  the  study  of  spinal  motions  must  be 
able  to  produce  simultaneous  motions  in  more  than  one  plane  and  to 
measure  loads  and  displacements  in  three  dimensions.  Liu,  et  al.  [18] 
subjected  cervical  FSU's  to  individual  modes  of  loading  (axial,  shear, 
bending,  and  torsion)  with  a  material  testing  machine  and  recorded 
motions  with  strain  gages.  The  Planar  Testing  Apparatus  used  here  is 
capable  of  producing  controlled  loading  or  three  components  of 
controlled  displacements  simultaneously.  The  PTA,  in  conjunction  with 
the  displacement  measuring  apparatus,  allows  for  three  dimensional 
motion  and  motion  measurement,  and  direct  data  analysis. 

The  three  dimensional  motion  measurement  system  and  analysis  of 
Eulerian  rotations  were  developed  and  refined  from  the  techniques 
presented  in  several  previous  studies  to  monitor  the  displacements  of 
the  FVS  specimens.  Tenser,  et  al.  [35]  used  a  displacement  measuring 
apparatus  very  similar  to  the  one  developed  here.  Measurement  errors 
due  to  the  assembly  design  and  limitations  on  the  range  capabilities 
of  his  transducers  were  reduced  in  the  improved  system  by  mounting 
transducers  in  a  gimbaled  arrangement. 

Methods  were  developed  here  for  the  precise  calculation  of 
translations  along  and  Eulerian  rotations  about  an  orthogonal  axis 
system  with  its  origin  at  the  center  of  the  moving  member  of  a  FVS. 
This  method  of  describing  three  dimensional  motion  of  one  body 
relative  to  another  is  simple  and  understandable.  Koogle,  et  al.  [14] 
used  the  kinematic  analysis  developed  by  Dieudonne,  et  al.  [6]  for 
the  calculation  of  Eulerian  rotations  of  his  test  fixture  but  had 
difficulty  with  the  elastic  strain  gages  used  to  monitor  displacements 
and  with  the  design  of  the  fixture  to  hold  spine  specimens.  The 
analysis  used  in  the  present  study  is  well  suited  to  our  displacement 
measuring  apparatus. 

Load-displacement  tests  of  seven  cervical  FVS's  yielded  several 
interesting  results.  There  was  no  dependence  of  stiffness  on 
displacement  rate  at  the  rates  used  in  this  study.  This  study 
suggests  that  velocities  greater  than  5  mm/sec  or  5  deg/sec  are 
required  to  observe  viscoelastic  effects.  Curves  showing  axial  and 


shear  forces  vs.  axial  and  shear  displacements,  and  moment  vs. 
rotation,  displayed  large  neutral  zones,  or  regions  of  small  load  and 
low  slope  as  displacements  increase.  All  specimens  were  stiffer  in 
compression  than  in  tension,  and  the  specimens  from  the  upper  cervical 
region  were  stiffer  in  response  to  shear  and  vertical  motions  than 
specimens  from  the  lower  region.  Specimens  taken  from  the  upper 
cervical  region  showed  greater  stiffness  to  posterior  motion  than  to 
anterior  motion. 


Few  previous  i_n 
performed,  and  the  t 
established.  For  th 
to  increasingly  high 
the  cervical  spine, 
bony  fracture,  and  i 
even  at  flexion  angl 
compressive  loads  of 
specimens.  It  appea 
structural  elements 
particular  structure 


vitro  investigations  on  cervical  spine  have  been 
olerance  of  this  region  to  injury  is  not  well 
is  reason,  cervical  FVS's  were  rotated  in  flexion 
er  angles  to  determine  the  fracture  behavior  of 
All  specimens  showed  soft  tissue  failure  before 
n  three  of  six  cases,  no  bony  fracture  was  evident 
es  of  30°.  In  addition  to  high  flexion  angles, 
150  N  or  more  were  necessary  for  fracture  of  the 
rs  that  there  was  a  gradual  giving  way  of  all 
,  instead  of  a  catastrophic  fracture  of  a 

,  for  the  experimental  condition  used. 


5.  SUMMARY  OF  PROJECT 


The  methodologies  developed  using  first  lumbar  and  then  cervical 
spine  segments  here  will  serve  as  a  basis  for  further  sagittal 


investigations  of  the  spine.  Tests  of  additional  cervical  spine 
specimens  will  establish  the  statistical  significance  of  the  results 
gathered  thus  far.  Future  tests  should  consider  more  severe  load 
combinations.  The  sagittal  results  presented  here  and  the  results 
from  these  future  tests  will  provide  a  comprehensive  description  of 
the  cervical  spine  and  its  tolerance  to  flexion  type  injuries.  The 
Planar  Testing  Apparatus,  displacement  measuring  device,  and  kinematic 
analysis  are  an  advanced  biomechanics  testing  system  suited  to  studies 
of  spine,  as  well  to  investigations  into  the  behavior  of  the  knee, 
wrist,  shoulder  or  other  complex  anatomic  joints. 
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Figure: 10  Test  stage  geometry 


V  :  Voltage  e  -  error 

c  -  coanand 
fb  -  feedback 

K  :  Aaplification  of  the  servo-controller 
L  :  Aaplification  of  the  LVDT 
P  :  Crossection  area  of  hydraulic  piston 
x  :  Piston  position 
a  :  Tiae  constant 
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11  Displacement  control  system  block  diagram 
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21.6  Creep  displacement  of  intervertebral  joint  at  three 
sagittal  rotational  angle  settings. 
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2.2  Creep  displacement  response  with  exponential  fit  of 
5-parameter  Kelvin  Solid  Model  in  parallel  of  the 
intact  intervertebral  joint.  (to:time  at  initial 
condi tion) 
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Fi gu re : 22 . 4  Creep  displacement  response  with  exponential  fit  of 
5-parameter  Kelvin  Solid  Model  in  parallel  of  the 
intact  intervertebral  joint,  (to: time  at  initial 
condi tion ) 


dlapl  (mm) 


\  ■%  *%r.  ^."^V.*V~. *V"  '■,"  *."  *\r  \r  v.”  *■.''  <*"  ^  VI 


kaiv&n  solid  aodal  at  zaro  aagittal 
angie  ISoecise”  5  to  -  19  sec) 


/ 


‘1 


‘  *r - — 1 - •- — - ' - 1 - 1 - - - 1 - 4 - 4— 

o  404  4"  90c  wo  idoc  uoe  i«e  teoo  woe 

Time  (sec) 

Figure: 22. 5  Creep  displacement  response  with  exponential  fit  of 
5-parameter  Kelvin  Solid  Model  in  parallel  of  the 
intact  intervertebral  joint.  (to:time  at  initial 
condi tion) 
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Load  supported  by  the  intervertebral  joint  at  three 
sagittal  rotational  angle  settings. 
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Fi gure : 24 . 2  Load  supported  by  the  i ntervertebral  joint  at  three 
sagittal  rotational  angle  settings. 
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ur'e:24.3  Load  supported  by  the  intervertebral  joint  at  three 
sagittal  rotational  angle  settings. 
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Figure:24.6  Load  supported  by  the  intervertebral  joint  at  three 
sagittal  rotational  angle  settings. 
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Figure: 29  Three  diaensional  coordinate  systea  for  the  Functional 
Vertebral  Segaent.  The  test  stage  coordinate  systea  has 
it's  origin  at  the  top  surface  and  center  of  the  stage 
with  axes  positive  in  the  saae  direction  as  the 
coordinates  shovn  in  the  figure. 


Figure:30  A  transducer  group  consisting  of  one  LVDT  and  two  RVDT 
and  supporting  structure 


31  Reference  axes  and  transducer  group  orientation  for 
aeasuraaent  of  the  displaceaent  of  the  aoving  vertebra 
vith  respect  to  tha  fixed  verebra.  O.aorigin  of 
transducer  group,  0M-origin  of  aovinglvertebra, 
O.aorigin  of  fixed  vertebra 


F i gu re : 32  Test  Fixture  (modified  from  Kou,  1935) 


Relationships  between  the  transducer  groups  origin,  A 
the  center  of  the  fixed  body,  S«,  the  center  of  the 
■oving  body,  S  ,  and  a  point  oncthe  specimen  collar,  A 


Z  Spec iaen  Displace 


Z  Spaciaen  Olaplacai 


DZ3ZBZVZ 


CZ2XZVZ 


39.2 


HZUXZVZ 


Specimen  Olaplaceaant 


eci 


CZ2XYZVZTS 


specimen  di  splacement  (x,y,z)  vs 
displacement  (z)  tor'  7  specimens 


Stag*  OliplactMnt 


HZUAXYZVTS 


BZ4XYZVTS 


Stags  Olsplacaaant 


Stags  Olaplacaaent 


iien 


Spaclaan  Oisplacai 


Diaplacai 


axaxisxvx 


Specimen 


ti-w 


/  /. 


eel 


CX3XZVX 


aceiiient  (a)  for  7  specimens 


hxsxzvx 


Oiaplaci 


EX2XZVX 


Specimen  Oisplace«ant 


*ec  i 


BX4XZVX 


OiipUcMant 


SpaciMn  Olapiacaaant 


Olapl 


X  Spaciaan  Oiaplacaaanfc 


Figures : 44 . 1 -44 . 7  Moment  (y)  vs  rotation  (theta)  for  7  specimens 
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Moment  (y)  vs  rotation  (theta)  at  large  angles  for 
3  specimens  -  2  different  angles  for  each  specimen 
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0 

Angle 

1  Max  vert 

2  Max  hori 

(  deg  ) 

range  (as) 

range  (am) 

1 

0 

127.00 

330. 20 

2 

1  0 

104.14 

274.32 

3 

20 

81.28 

223.52 

4 

30 

63.50 

167.64 

5 

40 

43.18 

114.30 

6 

50 

22.86 

58.42 

Table: 1  Range  of  motion  of  Planar  Testing  Apparatus 
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Angle 
(  deg  ) 
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2  Vert 
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3  Hoient 
(  N  -  •  ) 

1 
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28912-34694 

22240-28912 

1582 
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34250-36474 

18682-22685 

2712-2938 
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40 
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23  1  30-3  1  1  36 
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30 
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1  2  0  1  0-  2  3  1  30 
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0 
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Table :  2  Force  ranges  at  maximum  and  minimum  positions  for  each 
angle  along  the  vertical  axis  of  the  Planar  Testing 
Apparatus . 
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Table: 3  Force  ranges  at  maximum  and  minimum  positions  for  each 
angle  along  the  horizontal  axis  of  the  Planar  Testing 
Apparatus . 
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U  T  EDWARDS  ET  AL  1986  AFOSR-TR-86-0452 
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0  Test  ID  1  SLOPE 


1 

B 

0.002408 

2 

D 

0.006248 

3 

G 

0.002376 

4 

H 

0.001872 

2  INTERCEPT 


84.297449 

63.229373 

83.011050 

78.306532 


Table:6  Linear  regression  of  .  disc  load  at  zero  sagittal 


Sp«c. 

FVS 

ID 

Sex 

Weight 

Cause  of  death 

B 

C6-T1 

666638 

N 

66 

110 

COPD 

C 

C2-C4 

1647183 

M 

52 

163 

Cardioayopathy 

D 

C6-T1 

192343 

P 

76 

124 

CVA 

B 

C3-C5 

192543 

P 

76 

124 

CVA 

P 

C2-C4 

616886 

M 

73 

150 

CA  colon 

G 

C5-C7 

53852 

N 

69 

120 

Kyocard.  infarct 

H 

C2-C4 

53852 

N 

69 

120 

Hyocard.  infarct 

Table:7  Test  Specimen  Information 


APPENDIX  A 

Transducer  Calibration 

The  LVDT's  vere  calibrated  using  a  depth  gage  with  a  resolution 
of  0.25  mm.  The  transducer  group  housing  block  was  clamped  to  a 
benchtop  and  the  depth  gage  was  held  stationary  at  the  same  height  and 
coaxially  with  the  LVDT  support  rod.  The  gage  tip  was  placed  in  the 
cup  on  the  support  rods  end.  With  the  initial  location  of  the  rod  tip 
displaced  20  mm  from  its  position  corresponding  to  zero  voltage,  the 
gage  was  turned,  forcing  the  support  rod  back  towards  the  LVDT  body 
and  towards  the  zero  voltage  location.  Readings  were  recorded  from  a 
digital  voltmeter  every  2  mm  for  a  total  of  40  mm.  Voltage  vs. 
displacement  was  graphed  and  a  curve  was  fit  using  an  RS/1  software 
package  (RSI  software,  Bolt,  Beranek,  and  Newman  Research  Systems. 

The  RVDT's  were  calibrated  using  a  40°  sector  divided  by  radii  10  inch 
es  in  length  and  placed  2°  apart.  The  transducer  group  housing  block 
was  clamped  onto  the  benchtop  with  the  pivot  point  exactly  above  the 
origin  of  the  sector  and  the  LVDT  body  directed  at  90°  to  the  housing 
block.  The  LVDT's  support  rod  tip  was  aligned  with  each  radii  and  the 
voltages  recorded  via  a  digital  voltmeter.  Voltage  vs.  rotation  was 
graphed  and  a  curve  fit  using  RS/1.  The  calibration  curves  for  the 
LVDT's  and  the  RVDT's  are  given  in  A: 1-2. 


APPENDIX  B 


Kinematics 

The  components  of  the  equations  used  to  describe  the  kinematics 
are  described  in  detail  below. 
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where  C  stands  for  cosine,  S  stands  for  sine,  and  9  ...  is  the  angle  of 
rotation  of  the  transducer  groups  on  the  PTA  briage  about  the  Z^q 
axis.  This  is  referred  to  as  the  set-up  angle. 
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where  9y.  is  the  angle  of  rotation  of  the  transducer  group  pivot  point 
about  tne  Y  axis  as  determined  by  RVDT  2 
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where  8l.  is  the  angle  of  rotation  of  the  transducer  group  pivot  point 
about  tne  Z  axis  as  determined  by  RVDT  1 


where  D.  is  the  distance  (mm)  from  the  point  A  to  the  point  A  as 
determined  by  LVDT  1  0 
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where  ii,  9  and  <|>  are  successive  rotations  of  S  about  X  ,  Y  and  Z 
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where  fj(a)  is,  in  expanded  terms, 

,  ,  .  J  .2  2  „2  d2  „1  111 

f.(a)  =  A  +  A  i2  +  A  .j  +  B  .j  +  B  i2  +  fi  i3  +  x  +  y  +  z  - 

IMa  +  2(x  '  Bil^AilTll  +  Ai2T21  +  Ai3T31)  + 

2(y  -  Bi2)(AilTi2  +  Ai 2T2 2  +  Ai3T32)  + 

2(z  -  Bi3)(A.1T13  +  A.2T23  +  A.3T33) 


2(xB. .  +  yB.„  ♦  zB,7) 


vhere  implies  f^(a)  and  the  six  expanded  f^a^)  equations  are 


3f.(«)  =  2(x  +  A.1T11  +  A.2T21  +  A.3T31  Bn) 


3f.(a)  -  2(y  +  A.jT^  +  A.2T22  +  Ai3T32  -  B-2) 


af.(a)  =  2(z  +  A.1T13  +  A.2T23  +  A.3T33  B.3) 


3fj(ot)  =  -2(x-Bil)(AilT12  +  A.2T22  +  Ai3T32)  + 
3*  2(y  -  Bi2)(AilTll  +  Ai2T21  +  Ai3T31) 


3f.(a)  =  2(x  -  Bil)(-AilS0Ci|/  +  Ai2S«iCeC^  +  A.  3C<|>C0Cij/)  + 

30  2(y  -  Bi2)(-A.ises\j/  +  A.2S<|C0S^  +  A.3C<fCeSg/)  - 

2(Z  -  Bi3)(AilC9  +  A.2s<t>se  +  Ai3C4>S0) 

3fi*g)  =  2(x  "  Bil^Ai2T31  "  Ai3T21)  + 

3*  2(y  -  Bi2^Ai2T32  "  Ai3T22 ^  + 


2(z  B^a)(Ai?^n  A-ii^?'t) 


APPENDIX  C 


Error  Analysis 

To  determine  the  accuracy  of  the  equations  for  the  three 
translations  and  three  Eulerian  rotations,  along  with  the  accuracy  of 
the  displacement  hardware,  a  trial  displacement  collar  was 
constructed.  Three  balls  were  attached  to  a  lead  cylinder  which  was 
positioned  at  the  center  of  the  PTA  test  stage.  The  balls  were  at  a 
known  height  and  at  known  positions  with  respect  to  the  center  of  the 
stage.  The  cups  on  the  LVDT  support  rod  tips  were  linked  to  the 
balls,  exactly  as  if  the  collar  were  on  a  specimen  Figure  C : 1 .  The 
test  stage  was  displaced  by  three  mm  along  the  X  and  Z  axes  and  by 
three  degrees  about  the  Y  axis.  The  translation  and  rotation  data 
generated  by  the  displacement  software  were  compared  to  the  input 
data.  The  results  for  uniaxial  command  displacements  are  shown 
graphically  in  Figure  C:2.  Maximum  errors  of  0.4  mm  for  test  stage 
translations  of  3  mm  and  1.2  for  stage  rotations  of  10°  resulted.  The 
error  in  translation  affects  the  measured  rotation  angle. 
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